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CHAPTER  1 
INTRODUCTION 


The  ElectroScience  Laboratory  (ESL)  has  been  involved  in  the 
detection  of  underground  targets  of  one  type  or  another  for  over  a 
decade.  During  the  past  3  years,  part  of  our  efforts  have  been 
directed  towards  the  detection  of  tunnels  using  the  Video  Pulse  Radar 
concept  developed  and  used  previously  for  the  detection  of  such 
targets  as  utility  lines  and  geological  structures.  One  of  the  major 
problems  for  this  kind  of  radar  system  is  that  a  wide  variety  of 
unwanted  targets,  classified  here  as  clutter,  can  be  present  in  the 
received  waveform. 

A  major  effort  has  been  underway  at  Geo-Centers  to  provide  a 
target  identification  capability  through  mapping  techniques  using 
data  provided  by  ESL  from  Video  Pulse  Radar  Systems.  Due  to  the  need 
for  more  immediate  results  to  evaluate  the  progress  of  our 
measurement  program,  a  small  scale  effort  has  been  devoted  to  an 
in-house  mapping  project.  Obviously,  these  results  are  not  to  be 
considered  competitive  with  the  major  effort  but  since  they  contain 
interesting  findings,  they  are  the  subject  of  the  present  thesis. 

An  important  caveat  needs  to  be  included  at  this  point,  namely 
that  mapping  itself  may  not  be  possible  at  certain  sites  or  times 
because  of  restrictions  on  the  placement  of  the  antenna  on  the 
surface.  This  limitation  was  experienced  at  several  sites,  including 
mountainous  terrain  and  wooded  regions.  In  some  cases  only  a  single 
antenna  position  is  available. 

A  companion  study,  modifying  techniques  developed  under  contract 
DAAK70-77-C-0114  appears  to  have  a  capability  of  identifying  tunnels 
on  such  a  single-look  basis.  This  new  approach,  however,  cannot  be 
completed  under  the  present  contract  because  of  available  time  and 
funding;  but  some  preliminary  results  are  very  promising. 

Basically  the  radar  used  in  this  study  is  a  simple  video,  or 
base-band  pulse  system.  A  short  pulse,  of  either  150  picosecond  or  6 
nanosecond  duration,  is  connected  to  one  element  of  a  crossed,  folded, 
loaded  dipole  antenna  system.  This  pulse,  of  course,  has  a  very 
broad  spectrum,  extending  from  the  repetition  rate  of  the  pulser  to 
the  lower  microwave  region. 
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The  crossed  antenna  pair  acts  as  the  principal  filter  in  the 
system.  Its  response  has  been  broadened  by  loading  techniques  [1-3] 
and  may  be  modeled  in  the  complex  frequency  domain  by  a  single  pole 
pair  with  a  significant  real  component.  For  the  size  of  tunnel 
considered  iri  this  study,  the  crossed  dipole  antenna  will  pass 
frequencies  up  to  and  just  beyond  the  first  tunnel  resonance.  The 
crossed  antenna  configuration  also  eliminates  all  returns  from 
planar  parallel  surfaces.  Thus,  a  great  deal  of  processing  of  the 
signal  is  done  in  analog  fashion  by  the  radar  itself. 

Four  specific  antenna  systems  discussed  in  detail  by  Wald  L3] 
are  used  in  this  mapping  study.  They  are:  a)  a  Terrascan  antenna, 
b)  a  modified  Terrascan  antenna,  c)  a  short  box  antenna,  and  d)  a 
long  box  antenna. 

The  modification  to  the  Terrascan  antenna  consisted  of  moving 
the  balun  transformer  closer  to  the  antenna  terminals,  thereby 
shifting  multiple  reflections  from  it  out  of  the  desired  range 
window.  The  box  antenna  was  lengthened  just  prior  to  a  field  test 
of  the  antenna  because  of  computations  we  made  of  the  field  scattered 
by  a  square  tunnel.  These  computations  indicated  that  the  tunnel 
resonance  was  lower  in  frequency  than  we  had  expected. 

Data  was  obtained  for  mapping  studies  from  three  sources.  The 
first  of  those  sources  was  an  analytic  computation  made  by  Davis  [4] 
using  a  circular  cylinder  of  e  =1  immersed  in  a  homogeneous  medium  of 
er=9  and  a=0.01  as  a  model  of  a  tunnel.  The  other  two  sources  were 
field  trips  made  to  a  site  in  Gold  Hill,  Colorado  known  as  the 
Hazel  A  mine  and  also  to  the  Curtis  Schoolyard  in  Trumbull  County, 
Ohio  where  a  subsidence  problem  had  been  called  to  our  attention  by 
Professor  Anne  Harris  of  Youngstown  University.  Professor  Harris, 
a  geologist,  and  Mr.  Richard  Myers,  of  the  Bureau  of  Mines,  have 
both  been  very  helpful  to  us  in  our  studies. 

Before  developing  the  mapping  concept,  a  rather  naive  discussion 
of  the  scattered  time  domain  signal  from  a  tunnel  is  presented.  We 
model  the  tunnel  as  a  slab  with  e  =1  immersed  in  the  earth  medium  of 
f  >1  (see  Fig.  1).  A  pulse  is  rnormally  incident  on  the  slab  from 
tne  earth  and  results  in  a  reflection  from  the  upper  interface.  The 
corresponding  reflection  coefficient  is  positive.  A  pulse  is  also 
transmitted  into  the  slab  and  is  reflected  from  the  lower  interface. 
The  corresponding  reflection  coefficient  here  is  negative,  [he 
reflected  pulse  now  propagates  back  to  the  upper  interface  where  it 
is  both  transmitted  back  to  the  receiving  antenna  and  reflected  once 
again  into  the  slab.  Figure  1  illustrates  this  (note  that  the  rays 
are  actually  coll  inear  but  are  shown  separately  for  clarity).  These 
multiple  reflections  thus  set  up  a  train  of  reflected  pulses  that 
result  in  a  potential  multiple  lobe  signal  structure.  In  the  mapping 


process,  these  multiple  returns  will  show  up  as  i  train  of  apparent 
scatterers  at  depths  (from  Fig.  1J  related  to  2do,  2do+2h,  2do+4h, 
etc.  In  fact,  the  very  presence  of  such  a  train  in  a  map  is  an 
indication  that  a  target  whose  resonance  lies  within  the  bandwidth 
of  the  radar  system  is  present. 

As  has  been  indicated,  the  goal  of  this  thesis  is  to  develop  a 
mapping  technique  and  then  apply  it  to  experimental  data.  Since 
measurements  are  extremely  time  consuming,  there  are  a  limited  number 
of  positions  (as  stated  before)  at  which  an  antenna  may  be  placed. 

This  leads  to  a  necessary  compromise  between  the  number  of  data 
points  desired  for  the  construction  of  a  map  and  the  number  that  can 
actually  be  obtained  in  the  time  allotted  for  data  gathering. 
Obviously,  even  for  a  very  smooth  earth  one  cannot  afford  to  take 
data  with  antenna  positions  separated  by  one  centimeter  if  an  area  of 
say  one  square  kilometer  is  to  be  mapped.  Even  if  this  were  feasible, 
the  ground  conditions  would  very  likely  change  over  the  time  period  in 
which  the  data  was  being  aquired.  As  a  consequence,  the  mapping 
method  used  must  oe  able  to  deal  with  sparsely  sampled  areas. 

Fig.  2  illustrates  the  mapping  technique  to  be  used.  In  this 
illustration  we  have  6  waveforms  recorded  at  6  successive  locations 
on  the  surface.  These  are  displayed  with  the  time  axis  vertical 
in  Fig.  2.  Each  of  these  waveforms  is  stored  digitally,  consisting 
of  a  certain  number  of  samples.  The  format  of  our  final  grey  level 
plot  will  have  the  same  orientation  as  the  waveforms  shown,  namely 
the  horizontal  axis  will  be  traverse  along  the  ground  and  the 
vertical  axis  will  be  time  or  depth.  This  yields  a  cross-sectional 
view  of  the  subsurface.  We  will  be  generating  one  horizontal 
line  of  this  map  at  a  time. 

To  start  the  process,  we  pick  a  time  delay  (or  depth  if  the 
ground  parameters  are  known)  at  which  to  start  the  map.  This  is 
specified  as  the  waveform  sample  number  where  we  want  to  start.  This 
starting  sample  number  is  shown  as  the  horizontal  line  labeled 
"reference  depth"  in  Fig.  2.  One  sample  from  each  waveform  in  the 
traverse  (the  dot  on  the  reference  depth  line)  is  taken  at  this  time 
delay.  The  amplitude  of  these  samples  is  shown  plotted  again  on 
the  lower  axis  in  Fig.  2  (dots  on  lower  curve). 

We  could  at  this  point  make  a  map  from  just  these  waveform 
samples  but  it  would  be  extremely  small  and  hard  to  read.  Instead 
the  6  waveform  samples  are  fed  to  an  interpolator,  which  fits  a 
smooth  curve  through  them  as  shown.  We  may  interpolate  as  many 
points  as  we  wish  between  the  original  waveform  samples  to  yield  as 
large  a  map  as  we  want. 


Since  we  will  be  using  a  line  printer  as  our  output  medium  we 
do  not  have  a  continuum  of  grey  levels  available  to  represent  the 
amplitude  of  the  interpolated  signal.  For  our  current  installation 
the  limit  is  14.  Thus,  the  vertical  axis  shown  as  "grey  level"  In 
Fig.  2  must  be  quantized  into  14  equal  intervals.  After  this 
quantization,  the  interpolated  signal  is  scanned  from  left  to  right, 
and  graphic  characters  representing  the  grey  level  of  each 
interpolated  point  are  output  to  the  line  printer.  Thus,  one 
horizontal  line  of  the  map  is  printed. 

To  complete  the  entire  map,  we  return  to  the  original  digitized 
waveforms  (top  of  Fig.  2)  and  move  down  one  sampling  interval  on  each 
one.  This  results  in  a  new  reference  depth  line  just  one  sample 
below  the  old  one. 

The  interpolation  and  quantization  processes  are  repeated  and 
another  horizontal  line  is  output  to  the  line  printer  just  below  the 
previous  one.  We  repeat  this  procedure  of  moving  down  one  sampling 
interval  on  the  orginal  waveforms  until  we  have  mapped  as  much  of 
the  data  as  we  require. 

If  a  two-dimensional  grid  of  antenna  positions  was  used  to  record 
data,  another  type  of  grey  level  plot  may  be  produced.  By  sampling 
all  waveforms  in  the  grid  at  the  reference  depth,  a  "horizontal  slice" 
or  plan  view  may  be  made.  Both  the  plan  view  and  cross-sectional 
plot  are  discussed  in  detail  in  Chapter  2. 

It  can  be  seen  that  the  mapping  technique  presented  is  well 
suited  to  produce  maps  from  sparsely  sampled  grids.  The  inclusion  of 
interpolation  provides  a  map  that  is  large  and  smooth  enough  for  the 
eye  to  read. 

The  structure  of  the  remainder  of  this  thesis  will  be  as 
follows:  Chapter  2  presents  the  interpolation  function  used  in  the 
grey  level  plotting  package,  while  Chapter  3  discusses  the  grey  level 
plotting  routine  itself.  Chapter  4  presents  results  of  application 
of  the  grey  level  mapping  technique  to  analytic  and  experimental  data. 
Finally,  Chapter  5  discusses  conclusions  and  recommendations  for 
further  improvement  to  this  method  of  analysis. 


CHAPTER  II 
INTERPOLATION 


In  this  chapter  we  will  derive  the  interpolation  polynomial  used 
as  part  of  the  grey  level  plotting  package  and  discuss  the  effect  this 
method  of  interpolation  has  on  the  grey  le'el  plots  themselves. 

A  formal  statement  of  the  interpolation  problem  is: 


Given  a  continuous  f(x)  sampled  at  x=x  ,  x-j ,  •  •  • ,  x„  where  x^<xj 
for  i<j,  find  a  continuous  approximation,  F(x)  to  f(x)  in  the  interval 
Xo<x<Xn  such  that  for  ar.y  e>0 

I f (x)-F(x) |  £  e 

for  each  xc[xo,x  ].  (Note:  In  this  study,  x  denotes  the  position  of 
the  antenna  on  tne  surface  of  the  ground.  It  is  not  used  to  designate 
time  of  arrival  or  the  depth  of  the  target.) 

By  use  of  the  Weierstrauss  Theorem[5],  it  can  be  shown  that  F(x) 
exists  for  any  given  f(x)  and  e. 

After  finding  F(x),  the  interpolation  process  itself  consists  of 
finding  values  of  F(x)  at  points 

x1-<x<x^+]  i  =0 , 1 ,  •  •  •  ,n-l 

To  derive  the  approximation  polynomial,  we  define  Fn(x)  as  a  poly¬ 
nomial  of  order  n  and  write 

Fn(x)  =  l  akxk  (1) 

k=0 


where  n+1  is  the  number  of  sample  points. 

We  then  constrain  pn(x)  to  match  f(x)  at  all  of  the  sampled 
points  x=x0,xi , • • • ,xn  so  that 


Fn(xj)  =  f(xj )  j  =0, 1 , • • • ,n 


We  also  require  that  the  xj  be  distinct. 


(2) 


Using  (1)  and  (2): 
n  . 

I  ak  xi  =  f(xj)  J=0>1 , • • • »n 
k=0  K  J  J 

in  matrix  form: 


1 
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C\J  o 
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xn 
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f(x0) 
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X1 
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Yn 

X1 
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f(x1) 
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2 

x2 

x2 

a2 

f(x2) 

1 

xn 

*  2 
xn 

xn 

n_ 

_an_ 

_f(xn)_ 

—I 

Xa=f 


If  the  determinant  of  is  non-zero,  then  a  solution  exists  to 
Xa=f  and  is  unique.  The  determinant  DET  X^  is  well-known  and  is  the 
Vandermonde  Determinant.  Its  value  is  given  by 

i=n 

j=n-l 

Vvxl'"-'xn>  =  TT  (*1-Xj) 

i=l 
j=o 
i  >J 


Since  the  xk  are  distinct,  Vn  =  DET  X  f  0  and  a  unique  solution 
to  Xa=f  exists. 


It  can  be  shown  (in  a  development  presented  in  Appendix  A)  that 
the  general  form  of  (1)  is 


Fn(x)  1  x  x2 

f(x0)  1  xQ  x* 

f(x])  1  x1  x2 


(3) 


f  (xn ) 
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■*  * 


The  coefficients  a^  of  (1)  can  be  derived  by  expanding  the  above 
determinant  by  the  first  row. 


A  useful  form  of  (3)  results  when  we  expand  instead  by  the  first 
column,  yielding 


x  x 


C1  X1 


X1  X1 


Fn(x) 


-  f(x0)  . 


xn  xn 


x  x* 
ii  n 


+  t(xn> 


I'  Vl  <-1 


Fn(x)  =  l  f(xk)Ln(x) 
k=0 


where  L^(x)  is  the  kth  Lagrange  Coefficient  in  x  of  degree  n. 


Upon  examination,  the  following  properties  of  the  Lagrange 
Coefficients  appear: 


(1)  they  are  polynomials  of  degree  n  in  x 


(2)  for  j,k=0,l , • • • ,n  and  j^k 


LS<xj)'° 


(3)  for  j,k=0, n  and  j=k 
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since  we  have  constrained  Fn (xj )  =f (x j )  for  j=0,l,---,n. 

Let  us  now  proceed  to  derive  the  explicit  form  of  the  L^(x).  From 
properties  (1)  and  (2)  we  know  that  L{^(x)  is  an  nth  degree  polynomial 
with  roots  x=x0,xj , •  •  •  ,xk_-|  ,xk+] , ••  •  ,xn  so  that 


l£(x)  =  ak(x-x0)(x-x1)'--(x-xk_1)(x-xk+1)”*(x-xn) 
n 

«k  T T  (x-x ^ )  (5) 

K  j=0  J 

j^k 

where  ak  is  an  arbitrary  (real)  constant  to  be  determined. 

Using  property  (3)  from  above  we  find 

.  L‘(V  .  i 

ak  -  - - - -  -  — - - 

TT  (xk-Xj)  IT  (xk-xj) 

j=o  K  J  J=0 

j^k  j^k 


substituting  back  into  (5)  we  have 


Ln<*> 


TT  jz~ 

j=o  Uk 

j^k 


(x-x  ) 

h 


(6) 


which,  when  used  in  conjunction  with  (4)  results  in  our  formula  for 
interpolation. 

This  Lagrange  formula  for  interpolation  has  several  attractive 
features: 

(1)  the  L*^(x)  depend  only  on  x  and  the  constants  {x^}  (i.e., 
the  sampled  points  of  the  original  function  -  Note  that 
the  xj  may  be  arbitrary  (but  distinct)  values.  In  par¬ 
ticular,  they  need  not  be  equally  spaced), 

(2)  the  L^(x)  are  dimensionless, 

(3)  the  form  of  the  equation  for  the  L^(x)  is  invariant  under 
a  change  of  variable  x=a+hs  (i.e.,  invariant  under  trans¬ 
lation  by  "a"  and  scale  change  by  "h").  Stated  differently, 
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the  Lp(x)  depend  only  on  the  relative  positions  of  the 
points  involved  in  the  interpolation, 

(4)  the  original  sampled  points  are  preserved  and  are  not 
altered  as  in  other  forms  of  curve-fitting  (e.g.,  mean- 
square  approximations). 

Figure  3  shows  an  example  of  a  Lagrange  interpolation  where  the 
interpolated  function  was  evaluated  at  19  points  between  each  known 
point.  The  actual  waveform  is  shown  at  the  top  with  the  sampling 
points  (in  this  case  non-uniformly  spaced)  shown  by  "x's".  We  show 
also  a  straight-line  approximation  to  the  actual  waveform  since  the 
exact  behavior  of  the  actual  waveform  is  known  only  if  the  samples 
satisfy  the  Shannon  sampling  theorem  for  the  highest  significant 
frequency  component  present.  The  smoothness  of  the  interpolated  curve 
in  the  botton  of  Figure  3  is  apparent. 

However,  Figure  4  exhibits  the  one  observed  weakness  of  the 
Lagrange  (or  any  polynomial)  method.  Large  discontinuities  in  slope 
between  adjacent  points  (especially  near  the  ends  of  the  interval) 
cause  an  "overshoot"  or  "peaking"  phenomenon  which  can  cause  annoying 
false  targets  to  appear  in  a  grey  level  plot.  The  peaking  can  be  seen 
in  the  interpolated  curve  at  the  bottom  of  Figure  4  near  sampling 
points  9  and  10.  More  dense  sampling  would  be  the  best  way  to  remove 
this  characteristic. 

The  generalization  of  the  Lagrange  interpolation  formula  to  the 
two-dimensional  case  is  quite  simple.  As  shown  by  Steffensen[6],  it 
is  just  the  product  of  two  one-dimensional  interpolations,  so  that, 

Fmn(x’y)  =  f(xk’^)LmS(x«y) 

k=0  i- 0 

=  l  l  f(xk»yJLm(x)LS(y) 

k=0  1=0  £ 

After  plots  were  made  of  the  results  of  several  interpolations  of 
experimental  data  using  the  entire  data  set  in  the  interpolation  form¬ 
ula,  it  became  apparent  that  a  better  mathematical  model  for  including 
the  physical  properties  in  the  mapping  process  was  needed.  It  is 
intuitively  unsatisfactory  to  let  points  physically  far  removed  from  a 
section  of  the  map  influence  the  interpolation  of  that  section  (e.g., 
points  in  the  lower  right  of  a  map  should  have  little  significance  in 
the  interpolation  of  the  upper  left  portion  of  the  map).  For  this 
reason,  some  sort  of  "localized"  interpolation  was  sought.  Several 
methods  were  investigated. 


ACTUAL  WAVEFORM 


ORIGINAL  DATA 


Since  the  Lagrange  coefficients  depend  only  on  the  relative 
positions  of  the  points  involved  in  the  interpolation,  it  was  felt 
that  by  weighting  them  in  proportion  to  their  distance  from  the 
desired  interpolation  point  the  required  localized  characteristics 
could  be  obtained.  The  success  of  this  method  depends  on  the  unaltered 


Lagrange  coefficients  themselves  being  weighted  in  proportion  to  their 
distance  from  the  interpolation  point.  Figure  5  is  a  plot  of  the 

Lagrange  coefficients  for  a  ten-point  grid  (i.e.,  x0=l,  x-|=2, - , 

xg=10)  as  a  function  of  the  interpolation  point  x. 

The  coefficients  were  calculated  using  (6).  If  the  expected 
behavior  were  exhibited  in  these  graphs,  the  coefficients  nearest  to 
the  arrows  (the  position  of  the  interpolation  point,  x)  should  be  the 
largest,  with  magnitudes  dropping  off  as  we  move  away  from  the  arrow. 
Figure  5  shows  that  this  is  not  the  case.  The  Lagrange  coefficients 
are  not  a  clear  function  of  the  position  from  the  interpolation  point 
(see,  e.g.,  plot  for  x=1.75).  Thus,  a  localized  interpolation  method 
weighting  the  Lagrange  coefficients  is  not  easily  accomplished. 


(2)  Small  Interpolation  Grids 

With  this  method,  the  sampled  data  is  partitioned  into  small, 
adjacent  subsets.  Then,  each  subset  is  used  to  interpolate  within 
its  boundaries.  As  can  be  seen  from  Figure  6,  the  smoothness  of  the 
resulting  curve  is  very  dependent  on  the  method  used  to  partition  the 
data.  Tne  points  included  in  each  small  interpolation  are  shown 
enclosed  by  the  brackets  at  the  top  of  each  graph.  Each  interpola¬ 
tion  (encompassing  3  sample  points)  is  butted  against  the  next  in  the 
figure.  The  partitions  used  in  the  top  curve  produce  a  relatively 
smooth  waveform,  but  the  partitions  used  in  the  bottom  curve  produce 
large  slope  discontinuities  in  the  waveform.  These  discontinuities 
would  have  a  detrimental  effect  in  any  resulting  grey  level  plot. 


(3)  Overlapping  Small  Interpolation  Grids 

The  final  method  tested,  and  the  one  that  yielded  the  best 
results  was  similar  to  method  (2)  but,  where  different  partitions  over¬ 
lapped  they  were  averaged  together  to  yield  the  final  interpolated 
value.  The  result  of  averaging  the  overlapping  portions  of  Figure  6 
is  shown  in  Figure  7.  This  data  was  the  same  as  that  shown  in  the 
top  of  Figure  4  and  the  bottom  of  Figure  4  (using  a  full  data  set 
interpolation)  may  be  compared  to  Figure  7  ("localized"  interpolation). 
Such  a  comparison  shows  that  the  localized  interpolation  gives  a  curve 
with  much  less  "overshoot"  -  particularly  noticable  around  points  2,  4, 
and  9,  and  yet  is  still  smooth.  This  result  is  to  be  expected  since 
lower-order  polynomials  are  used  in  the  localized  method  and  very  fast 
slope  changes  cannot  occur. 


CHAPTER  III 

GREY  LEVEL  PLOTTING  ROUTINE 
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The  primary  function  of  the  grey  level  plotting  routine  (here¬ 
after  GLPR)  is  to  take  an  array  of  data,  interpolate  it  to  expand  it 
to  a  given  size,  quantize  the  resulting  array  into  a  specified  number 
of  levels,  and  then  plot  that  information  in  the  form  of  grey  levels. 

Our  use  of  GLPR  has  been  primarily  to  plot  groups  of  experi¬ 
mental  underground  radar  returns.  The  resulting  grey  level  plot  is 
then  used  as  an  indication  of  the  radar's  ability  to  detect  a  target 
of  certain  prescribed  characteristics. 

We  will  now  discuss  the  various  inputs  to  GLPR.  The  array  of 
known  data  (the  "input  array")  may  be  constructed  in  several  ways, 
depending  on  the  method  used  to  gather  the  data.  If  a  two-dimensional 
grid  of  data  is  available,  this  may  be  used  directly  as  the  input 
array.  Subsequent  interpolation  and  plotting  will  yield  "plan  views" 
of  the  area  mapped.  Since  using  only  one  point  from  each  waveform  at 
the  depth  of  interest  for  the  entry  in  the  input  array  would  be  mis¬ 
leading  (because  the  equipment  used  in  data  taking  may  not  have  that 
high  a  frequency  response)  another  method  for  calculating  the  entries 
in  the  input  array  is  used.  A  depth  window  (i.e.,  time  window  on  the 
waveform)  is  set  up  on  each  waveform  and  the  RMS  value  of  the  curve 
calculated  therein.  (Choice  of  the  window  size  is  critical.  For  a 
discussion  of  this  point  see  Chapter  4.)  The  input  array  then  con¬ 
sists  of  these  RMS  values  for  a  given  depth  of  interest. 

If  only  a  one-dimensional,  or  linear  scan  is  available,  the 
input  array  may  be  set  up  using  all  sample  points  from  each  waveform 
as  a  column  in  the  array.  This  format  is  similar  to  that  used  in 
seismograms  [7].  In  order  to  keep  the  plots  to  manageable  size,  only 
a  portion  of  the  waveforms  in  the  time  window  of  interest  is  usually 
plotted. 

A  summary  of  all  array  conventions  used  in  this  section  is  shown 
in  Figure  8. 

It  is  not  necessary  for  the  map  data  to  be  equally  spaced.  GLPR 
will  allow  non-uniform  spacing  between  known  points  and  will  inter¬ 
polate  to  correct  for  it.  Although  the  spacing  may  be  non-uniform, 
it  must  be  consistent  throughout  the  map.  To  understand  the  differ¬ 
ence  between  consistent  and  inconsistent  non-uniform  maps,  see  Figure  9. 
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ARRAY  CONVENTIONS 
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O  =  KNOWN  DATA  POINT 


SUBSCRIPT  FORM  (Y,X) 
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VERTICAL:  Y -DIRECTION 
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O  O  O 
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3X3  ARRAY 


ARRAY  SIZE:  ROW  X  COLUMN 


Figure  8.  Statement  of  Array  Conventions. 
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INCONSISTENT  NON-UNIFORM  MAP 
(  POINTS  NOT  LINED  UP  ) 

Figure  9.  Consistent  and  Inconsistent  Non-Uniform  Maps, 
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Figure  10  shows  two  examples  of  physical  spacing  versus  array 
representations.  On  the  left  is  a  representation  of  the  physical 
layout  of  a  two-dimensional  map.  In  the  uniform  case,  spacing  be¬ 
tween  all  samples  is  equal  and  the  input  array  subscripts  reflect 
the  physical  position  of  the  sample  points.  In  the  non-uniform  case, 
however,  spacing  between  samples  in  the  horizontal  dimension  varies 
as  a  function  of  x  and  the  array  subscript  no  longer  reflects  the 
true  physical  spacing  of  the  sample  points.  Another  array  with  this 
spacing  information  must  be  supplied  to  GLPR. 

At  the  right  of  Figure  10  are  the  array  representations  of  the 
input  data  after  interpolation  (throughout  this  chapter  a  white  circle 
will  be  used  to  represent  sample  (known)  points  and  a  black  circle 
will  designate  an  interpolated  point).  The  size  of  the  interpolated 
array  is  specified  by  the  user  o'  GLPR  in  the  number  of  interpolated 
points  desired  between  each  known  point  (i.e.,  between  each  entry  of 
the  input  array).  At  the  top  right  of  Figure  10,  the  user  has 
specified  two  points  to  be  interpolated  between  known  values  in  the 
x-direction  and  one  point  between  values  in  the  y-direction. 

In  the  non-uniform  case  the  interpolation  is  handled  differently. 
The  interpolation  attempts  to  fill  in  "missing"  data  points  to  make 
the  final  plot  look  true  (i.e.,  without  the  distortion  that  might 
result  if  point  (1,2)  in  Figure  10  was  considered  to  be  equidistant 
from  points  (1,1)  and  (1,3).  The  result  is  that  interpolation  beyond 
that  specified  by  the  user  is  performed.  In  the  lower  right  of 
Figure  10,  the  user  has  specified  one  point  to  be  interpolated  between 
known  values  in  the  y-direction  and  no  points  between  values  in  the 
x-direction  (as  seen  from  the  lack  of  an  interpolated  point  between 
columns  1  and  2).  Note  that  interpolation  was  still  performed  in  the 
x-direction,  however,  to  make  the  final  plot  reflect  the  physical  lay¬ 
out  of  the  data.  Note  also  that  the  input  array  may  be  non-uniform 
in  one  direction  only  or  in  both. 

In  order  to  fully  explain  the  algorithm  used  by  GLPR,  several 
terms  must  now  be  defined. 

Block:  A  block  is  the  area  bounded  by  four  adjacent  points  of 
known  data.  Figure  11  shows  two  examples  of  blocks:  niform  on  the 
left  and  non-uniform  on  the  right. 

Grid:  The  pattern  of  measurements  (experimental  data)  that 
constitute  a  map. 

Subgrid:  A  subgrid  is  the  rectangular  set  of  known  points  used 
to  interpolate  inside  its  boundaries.  It  has  dimensions  VxX.  A  sub¬ 
grid  is  identified  by  its  upper  left  noint.  Figure  12  shows  a  4x5 
input  array  with  a  3x3  subgrid  consisting  of  the  blocks  numbered  1 
through  4. 
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PHYSICAL  CONFIGURATION 
(WITH  INPUT  ARRAY  SUBSCRIPT) 
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Figure  10.  Expansion  of  Input  Array  by  Interpolation. 
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Overlap:  An  overlap  occurs  if  a  fixed  point  (Y,X)  is  contained 
in  more  than  one  subgrid..'  In  Figure  13,  if  we  are  using  a  3x3  sub¬ 
grid,  then  each  point  in  the  block  shown  will  be  overlapped  by  the 
subgrids  marked  with  Overlaps  other  than  those  indicated  will 

also  occur.  Note  that  the  line  AB  is  also  overlapped  by  the  subgrid 
at  A  and  that  the  point  B  is  overlapped  yet  again  by  the  subgrid 
identified  by  point  B.  These  types  of  overlap  are  very  cumbersome 
to  deal  with  computationally  and  are  ignored  by  GLPR.  The  conse¬ 
quences  of  this  are  discussed  below. 

"Interpolating  a  Block  with  Subgrid  (Y,X)"  means  to  interpo¬ 
late  values  at  all  points  in  a  block  using  the  subgrid  (Y,X). 

"Filling  a  Block"  means  to  calculate  a  value  for  each  point 
inside  a  block  by  averaging  the  overlapping  interpolated  values  at 
that  point. 

A  "Horizontal  Scan"  is  the  filling  of  all  blocks  from  left  to 
right  for  some  fixed  Y.  In  Figure  14,  for  the  3x5  input  array  shown, 
filling  blocks  1-4  constitutes  a  horizontal  scan. 

The  algorithm  used  in  GLPR  is  then  to  perform  horizontal  scans 
from  top  to  bottom  until  all  blocks  in  the  input  array  have  been 
filled.  This  methodology  means  that  at  points  where  "higher-order" 
overlaps  occur  (see  "overlap"  above),  the  interpolations  using  points 
to  the  right  and  downward  are  preferred  and  in  fact  overwrite  the 
results  from  other  interpolations  on  the  left  and  upward.  The 
detrimental  effects,  if  any,  of  this  method  have  not  yet  been  observed 
in  practice. 

We  now  state  a  procedure  for  parameter  choice  for  the  user  of 
GLPR.  After  statement  of  the  procedure,  two  examples  of  its  use  are 

')  ■»*  /pr 

I)  Choose  the  dimensions  of  the  input  array.  One  of  the 
dime  sions  must  be  less  than  or  equal  to  126,  as  this  is  the  width 
of  a  line  printer  page.  If  a  two-dimensional  map  of  data  was  taken, 
the  dimensions  of  the  input  array  are  just  the  dimensions  of  the  area 
of  interest  in  the  map.  If  a  linear  scan  was  taken,  the  usual  choice 
for  the  y-dimension  is  the  number  of  sample  points  in  the  waveforms 
and  for  the  x-dimension  is  the  number  of  waveforms  in  the  scan. 

II)  Choose  the  number  of  points  to  be  interpolated  between  the 
given  points.  This  number  is  specified  independently  for  the  x-  and 
y-dimensions.  When  choosing  these  parameters,  it  must  be  remembered 
that  in  the  non-uniform  case,  even  if  no  points  are  desired  between 
given  points,  the  GLPR  will  still  interpolate  enough  points  between 
the  given  ones  to  eliminate  any  distortion  that  might  otherwise  result. 

As  an  example,  suppose  waveforms  were  recorded  at  positions  1, 

5,  6  and  8  in  some  reference  grid.  GLPR  would  interpolate  values  for 
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positions  2,  3,  4  and  7  on  the  grid  even  if  the  user  specified  no 
interpolated  points  were  to  be  generated  between  known  values.  Con¬ 
versely,  if  the  user  had  specified  1  point  to  be  generated  between 
known  values,  the  program  would  then  calculate  values  for  positions 
1.5,  2,  2.5,  3,  3.5,  4,  4.5,  5.5,  6.5,  7  and  7.5  on  the  reference 
grid. 

III)  By  completing  step  II,  the  dimensions  of  the  interpolated 
array  have  now  been  fixed.  They  may  be  calculated  from  the  following 
formulas: 

For  the  uniform  case: 

interpolated  array  x-dimension  = 

total  number  of  points  generated  in  x-direction  = 

(x-dimension  of  input  array  -1) 

x  (number  of  interpolated  points  desired  in 

x-direction  +1)  +1 

interpolated  array  y-dimension  = 
number  of  points  generated  in  y-direction 
in  one  horizontal  scan  = 

number  of  interpolated  points  desired  in  y-direction  +2 

For  the  non-uniform  case: 

Calculate  x-size  =  highest  coordinate  in  reference  grid 
-  lowest  coordinate  in  reference  grid 

then  interpolated  array  x-dimension  = 
total  number  of  points  calculated  in  x-direction  = 
x-size  x  (number  of  interpolated  points  desired  in 
x-direction  +1)  +1 

Calculate  y-size  =  maximum  spacing  of  any  two 
adjacent  points  on  the  reference  grid. 

then  interpolated  array  y-dimension  = 

maximum  number  of  points  generated  in  y-direction 

in  any  one  horizontal  scan  = 

y-size  x  (number  of  interpolated  points  desired 

in  y-direction  +1)  +1 

It  will  be  noted  that  the  entire  interpolated  map  is  never  kept 
in  memory  at  one  time.  Only  the  data  from  one  horizontal  scan  resides 
in  core  at  any  time.  This  facilitates  processing  of  large  maps. 

IV)  Choose  the  number  of  grey  levels.  A  maximum  of  fourteen 
levels  may  be  specified.  The  representation  of  these  levels  on  a  line 
printer  page  is  shown  in  Figure  15.  Maximum  resolution  is  achieved 
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Figure  15.  Line  Printer  Representation  of  Grey  Levels. 
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with  a  large  number  of  grey  levels  and  maximum  contrast  is  obtained 
with  a  small  number,  so  a  compromise  is  usually  necessary.  (Also  see 
VI  below  "clipping"). 

V.  Choose  the  subgrid  size  to  be  used  for  interpolation.  The 
subgrid  size  in  the  x-  and  y-dimensions  are  specified  separately.  If 
no  interpolation  is  desired  in  a  given  direction,  that  subgrid  size 
is  specified  as  1  (and  the  number  of  desired  interpolated  points  in 
that  direction  (step  II)  should  be  declared  zero). 

The  size  of  the  subgrid  may  be  chosen  up  to  the  size  of  the  in¬ 
put  array  itself.  Usually,  to  model  the  physical  processes  involved, 
the  entire  data  set  is  not  used. 

In  the  non-uniform  case,  data  physically  far  removed  from  cer¬ 
tain  points  may  unduly  influence  the  interpolation  around  that  point, 
so  care  must  be  taken  when  choosing  the  subgrid  size. 

VI.  Specify  the  clipping  level.  One  method  of  increasing 
resolution  in  a  given  area  of  a  grey  level  plot  while  sacrificing 
detail  in  others  is  by  "clipping"  the  data.  That  is,  data  in  the 
interpolated  array  are  tested;  if  a  data  value  is  greater  than  or 
equal  to  a  specified  threshold  value  (the  clip  value)  it  is  set  equal 
to  the  threshold  value.  This  will  sacrifice  detail  in  the  clipped 
areas  and  will  allot  more  grey  levels  to  the  remainder  of  the  map. 

In  practice,  one  plot  is  usually  made  with  an  arbitrarily  high  clip 
value,  in  order  to  determine  what  the  most  effective  clip  value  will 
be. 

In  order  to  clarify  the  use  of  the  above  procedure,  the  follow¬ 
ing  two  examples  are  presented: 


EXAMPLE  1 

We  have  17  waveforms  of  256  samples  each  which  were  taken  along 
a  straight  line  at  intervals  of  1  foot.  The  time  window  seen  by  each 
waveform  is  200  ns  and  the  ground  in  which  the  measurements  were 
taken  had  a  relative  permittivity  of  approximately  9. 

Since  this  is  a  linear  scan,  we  decide  to  do  a  cross-sectional 
view  grey-level  plot,  using  all  17  waveforms  and  the  entire  200  ns 
time  window.  We  will  plot  depth  down  the  page  and  horizontal  position 
across  the  page.  The  dimensions  of  the  input  array  will  then  be 
Y=256,  X=17.  Since  the  data  was  taken  on  a  uniform  grid,  we  need  not 
specify  the  spacing. 

Suppose  that  we  wish  to  plot  the  data  letting  the  width  and 
height  of  one  grey  level  character  represent  the  same  physical 
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reference  dimension.  That  is,  we  wish  to  make  the  horizontal  and 
vertical  scales  the  same.  This  will  result  in  a  small  amount  of 
visual  distortion  (since  the  actual  width  and  height  of  a  line  printer 
character  are  not  equal)  but  is  a  useful  first  approximation. 

The  physical  dimensions  of  the  area  represented  by  our  grey 
level  plot  are: 

horizontal:  (17  waveforms  -1)  *1  foot/waveform  *0.305 
meter/foot  =  4.88  meters 

vertical:  (20P-n.O-l0-^|  rce ter/ns j_  =  20  meters 

but  this  is  the  round-trip  distance,  so 
the  actual  range  is  20/2  =  10  meters 


Let  us  let  one  character  on  the  grey  level  plot  represent  one 
sample  from  the  waveform  in  the  y-direction  (i.e.,  no  interpolated 
points  in  the  y-direction) Then  the  unit  of  length  represented  by 
one  character's  height  or  width  will  be 

unit  length  =  ^pofSs  =  °'°391  meters/P°int- 
Then  we  should  have 


4.88  meters 


0.0391  meters/point 


=  125  points 


on  the  horizontal.  We  can  come  close  to  this  number  by  specifying 
that  we  want  6  points  interpolated  between  each  known  value  in  the 
horizontal  direction.  We  would  then  have  a  plot  (see  formula  in 
Step  III) 

(17-1 )*(6+l )+l  =  113  characters  wide 

This  is  a  distortion  of  our  original  plan  for  the  plot  by  compression 
of  about  9.6%.  We  will  consider  that  to  be  acceptable. 

We  have  now  specified  the  desired  number  of  interpolated  points, 
namely,  Y=0,  X=6  and  the  interpolated  array  dimensions  are  (see 
formula  in  Step  III):  Y=0+2=2  and  X=1 13. 
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We  will  initially  get  as  much  resolution  as  possible  by  specify¬ 
ing  fourteen  grey  levels  and  an  arbitrarily  high  clip  value. 

In  choosing  a  subgrid  size,  we  will  choose  to  interpolate  in  the 
x-direction  only,  since  we  have  a  fairly  high  number  of  samples  in 
the  y-direction.  We  do  this  by  specifying  the  y-dimension  subgrid 
size  as  one.  We  choose  the  X  subgrid  size  to  be  3,  as  this  is  the 
smallest  value  that  still  allows  overlapping  to  occur. 


EXAMPLE  2 

We  have  a  grid  of  waveforms  available  in  the  configuration  shown 
in  Figure  16.  In  one  direction  the  measurements  are  spaced  at  half¬ 
foot  intervals  and  in  the  other  direction,  measurements  were  taken  at 
1,  3,  and  2  foot  intervals. 

We  choose  to  display  this  information  in  a  series  of  plan  views 
at  selected  depth  windows.  Let  us  assume  that  some  window  of  interest 
has  been  chosen  and  RMS  values  of  the  waveforms  in  this  window  are 
available. 

After  interpolation,  we  expect  that  the  non-uniform  dimension 
will  be  the  largest,  and  therefore  choose  to  plot  this  down  the  page 
(y-direction).  The  dimensions  of  the  input  array  are  then  X=5,  Y=4. 

Since  the  data  are  non-uniformly  sampled,  we  must  specify  the 
spacing.  In  the  x-direction,  the  samples  are  uniformly  spaced,  so 
that  we  may  say  the  measurements  were  taken  at  positions  x=l,2,3,4,5. 
In  the  y-dimension,  we  choose  1  foot  as  our  basic  unit  and  specify 
that  the  measurements  were  taken  at  positions  y=l ,2,5,7. 

Let  us  once  again  attempt  to  present  the  data  in  "true"  form  by 
making  the  x-  and  y-scales  on  the  plot  equal.  In  the  x-direction, 
specify  14  points  to  be  generated  between  known  values.  (This  was 
chosen  solely  as  an  attempt  to  fill  the  width  of  the  page.)  The 
formula  (Step  III)  for  the  interpolated  array  x-dimension  then  gives 

(5-l)*(14+l)+l  =  61  points 

Let  us  also  try  to  correct  for  the  fact  that  the  height/width 
ratio  of  a  line  printer  character  is  5/3.  Expanding  the  x-scale  by 
5/3  would  give  101.67  points.  We  can  come  close  to  this  by  changing 
our  desired  number  of  interpolated  points  to  24  instead  of  14.  Then 
the  interpolated  array  x-dimension  is  (5-1 )*(24+l )+l  =  101  points. 

To  make  the  interpolated  data  "true."  we  will  need  twice  as 
many  interpolated  points  in  the  y-direction  (twice  the  number  of 
points  needed  before  5/3  correction).  This  is  because  our  unit  of 
measurement  on  the  y-axis  was  twice  that  used  in  the  x-axis  (i.e.,  1 
ft.  vs.  0.5  ft).  Using  the  formula  from  Step  III  for  the  non-uniform 
case,  we  have: 
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Figure  16.  Layout  of  Measurement  Positions  Used  in  Example  2. 


Y-size  =  MAX{ 2-1 ,5-2,7-5} 

=  MAX{ 1 ,3,2} 

=  3 

interpolated  array  y-dimension  = 

3*(14*2+1)+1  =  88  points 

We  then  choose  to  obtain  14  grey  levels;  we  set  clipping  at  an 
arbitrarily  high  value  and  choose  a  subgrid  with  Y=3,  X=3  for  inter¬ 
polation  to  obtain  overlap  in  both  dimensions. 


CHAPTER  IV 
GREY  LEVEL  PLOTS 


In  this  section  we  will  present  results  of  applying  the 
previously  discussed  grey  level  plotting  program  to  experimental 
underground  radar  returns.  We  will  discuss  separately  data  taken 
at  several  locations  with  different  antenna  types  and  different 
pulsers. 

The  first  location  from  which  data  was  processed  was  the 
Curtis  School  in  Trumbull  County,  Ohio.  The  site  was  a  schoolyard 
with  an  abandoned  coal  mine  structure  underneath.  Subsidence  had 
already  occurred  in  areas  of  the  school's  playground.  The  depth 
of  the  mine  structures  in  this  yard  was  approximately  2.4  m.  The 
yard  has  subsequently  been  drilled  using  a  square  grid  with  7.62  m 
spacing  between  holes.  The  rock  composition  was  found  to  be  rather 
complex  and  included  sandstone,  shale,  and  much  broken  rock. 

Several  visits  have  been  made  to  the  site.  The  data  presented 
here  is  from  the  first  trip,  made  in  January,  1978.  The  temperature 
had  been  below  freezing  over  an  extended  period  and  a  deep  snow 
had  fallen.  However,  the  schoolyard  was  plowed  by  local  officials 
prior  to  our  measurements. 

Initially,  the  area  was  scanned  using  the  Terrascan  system. 
Based  on  signal  returns  observed,  a  portion  of  the  yard  was  selected 
for  detailed  study.  Data  was  taken  according  to  the  plan  shown 
in  Figures  17  and  18.  Figure  17  shows  the  general  layout  of  the 
yard  with  an  indication  of  the  area  mapped.  Figure  18  gives  the 
coordinates  used  to  identify  each  waveform  recorded. 

The  measurements  were  made  using  an  unmodified  Terrascan 
crossed-dipole  antenna  with  a  150  ps  pulser.  The  receiver  was 
a  Tektronics  sampling  oscilloscope.  Data  was  recorded  using  a 
microprocessor  system  (_8j  and  floppy  disk  recorder.  Waveforms 
obtained  using  this  system  are  shown  in  Figure  19.  Each  waveform 
essentially  shows  the  received  signal  level  as  a  function  of  time. 

The  first  type  of  plot  attempted  was  a  plan  view  of  the 
area.  The  waveforms  were  examined  and  a  time  window  from  approxi¬ 
mately  30  to  40  ns  was  chosen  as  the  depth  at  which  the  largest 
returns  were  occurring. 
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Numbers  derived  from  the  peak-to-peak  signal  amplitude  were 
calculated  in  this  window  for  each  position  shown  in  Figure  18 
and  the  values  are  shown  in  Figure  20.  Only  the  waveforms  from 
columns  5-17,  odd  and  rows  A,B,D,F,H,K  were  used  in  the  grey  level 
plot.  This  is  the  largest  set  of  positions  that  constitutes  a 
consistent  map  (see  Chapter  III). 

It  should  be  noted  that  at  the  time  this  plot  was  made,  the 
GLPR  could  not  handle  non-uniform  grids  of  data.  Thus,  the  decision 
to  include  rows  A  and  K  in  the  plot  results  in  inherent  distortion. 
GLPR  treats  rows  A  and  B  and  also  rows  H  and  K  as  equidistant  from 

each  other.  These  rows  are  also  plotted  in  the  same  manner  so 

that  the  resulting  plot  is  compressed  at  the  top  and  expanded  at 

the  bottom.  Also,  this  early  version  of  GLPR  could  only  perform 

interpolation  on  entire  data  sets,  hence  the  interpolation  is  not 
localized. 

The  grey  level  plot  itself  is  shown  in  Figure  21.  Contour 
lines  were  added  to  help  distinguish  between  the  different  graphic 
characters  and  to  observe  the  characteristics  of  the  interpolation. 

Next,  computer  calculation  of  RMS  values  in  a  time  window 
of  specified  position  and  width  was  added  to  the  GLPR.  This  relieved 
the  burden  of  hand  calculation  of  signal  amplitudes.  GLPR  was 
also  upgraded  to  handle  localized  interpolations  (as  discussed 
in  Chapter  II). 

Using  the  same  waveforms  as  above  (i.e.,  A,B,D,F,H,K,  5-17  odd) 
plan  views  were  plotted  using  a  50-point  time  window  for  calculation 
of  the  RMS  value.  Each  waveform  consisted  of  256  "points"  or  samples 
and  spanned  100  ns  in  time.  This  gives  a  scale  of  (100  ns)/(250 
pts.)  =  0.391  ns/pt.  So  a  50-point  window  would  be  approximately 
20  ns  wide.  This  window  was  scanned  across  the  waveform  set  at 
25  point  (vrIO  ns)  increments.  The  results  are  shown  in  Figures 
22-30. 

It  should  be  remembered  in  viewing  these  plots  that  only 
42  numbers  were  input  to  the  interpolator,  namely  the  RMS  values 
calculated  for  the  waveforms  recorded  at  the  intersection  of 
coordinate  lines  A,B,D,F,H,K  and  5-17  odd.  Only  these  42  points 
on  the  grey  level  plot  are  experimental.  All  other  points  were 
generated  by  the  interpolation  in  trying  to  fit  a  smooth  2-dimensional 
surface  of  the  specified  size  to  the  given  data  points. 

An  interpolation  grid  size  (see  Chapter  III)  of  3x3  was  chosen 
for  these  plots  and  this  resulted  in  averaging  of  at  most  4  overlaps 
at  each  point  in  the  plot.  It  sould  be  mentioned  that  the  compres¬ 
sion/expansion  distortion  mentioned  above  is  still  present  in 
these  plots  as  the  GLPR  had  not  yet  been  modified  to  deal  with 
non-uniform  data  sets. 
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Figure  22.  Curtis  School  Site  Plan  View. 
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Figure  26.  Curtis  School  Site  Plan  View. 
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Looking  at  the  plots,  it  is  apparent  that  the  area  of  interest 
in  all  (except  plots  for  0-19.5  ns  and  10.2-29.3  ns)  is  the  upper 
right  corner.  This  corresponds  to  the  area  around  position  A15 
in  Figure  18.  Of  particular  note  is  the  plot  for  68.8-87.9  ns 
where  quite  linear  edges  are  present,  indicating  the  possible 
presence  of  a  tunnel  at  this  depth. 

A  rough  estimate  of  the  ground  permittivity  in  the  area  is 
6  £  e  _<  9.  Using  e  =  7.5  we  find  that  a  delay  of  69  ns  corre¬ 
spond^  to  a  depth  ofr3.7  m  and  a  delay  of  88  ns  to  a  depth  of 
4.9  m.  This  agrees  with  the  distance  to  the  floor  of  the  tunnel 
of  3.7  m  as  determined  by  drilling.  The  top  of  the  tunnel  was 
at  2.4  m,  corresponding  to  44  ns.  Examining  the  grey  plots  containing 
this  depth  (29.7-48.8  ns  and  38.5-58.6  ns)  we  note  strong  signal 
returns  in  the  A15  position,  but  these  are  visible  in  other  plots 
also.  It  appears  that  the  tunnel  response  is  contained  in  this 
data.  However,  many  other  features  are  also  present.  We  believe 
that  this  clutter  is  caused  by  inhomogenities  in  the  tunnel  over¬ 
burden  due  to  changes  in  water  flow  during  the  nearly  100  years 
the  tunnel  has  been  in  existence.  Figure  19  shows  the  large  signal 
returns  that  occur  at  these  shallow  depths.  It  was  these  large 
returns  that  initially  motivated  mapping  of  this  area. 

Figure  31  shows  the  grey  level  distributions  of  the  plots 
of  Figures  22-30.  In  each  plot,  the  abscissa  is  the  number  of 
the  grey  level  from  1  (white)  to  14  (black)  and  the  ordinate  is 
the  percentage  of  the  total  map  occupied  by  that  grey  level.  It 
is  interesting  to  note  that  the  grey  level  plot  that  showed  the 
possible  presence  of  a  tunnel  most  clearly  (68.8-87.9  ns)  is  also 
that  which  has  a  more  nearly  uniform  grey  level  distribution. 

Perhaps  the  most  critical  parameter  in  processing  the  plan 
view  type  of  grey  level  plot  is  the  choice  of  window  size  for  RMS 
value  calculation.  If  a  window  is  chosen  too  large,  a  target 
(which  usually  appears  as  a  ringing  in  the  waveform)  could  acci¬ 
dentally  be  removed.  On  the  other  hand,  a  window  size  too  small 
may  represent  a  time  resolution  in  excess  of  the  frequency  bandwidth 
of  the  radar  system  used  to  take  the  data. 

Another  form  of  grey  level  plot  is  the  cross-sectional  view. 

The  cross-section  differs  from  the  plan  view  in  that  it  is  not 
usually  necessary  (or  even  desirable)  to  use  the  RMS  windowing 
process  discussed  above.  In  this  type  of  plot,  the  vertical  axis 
consists  of  a  section  of  the  input  waveform  plotted  point-by-point. 

An  absolute  value  of  the  waveform  is  taken  so  that  both  positive 
and  negative  peaks  in  the  signal  will  appear  on  the  grey  level 
plot  as  darkened  areas. 
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In  the  cross-section  plot  the  "signature",  or  shape  that 
a  target  assumes  in  the  plot  is  not  an  actual  cross-sectional  view 
of  the  target  itself.  As  shown  in  Figure  32,  at  each  point  in 
a  linear  traverse  the  distance  to  the  target  (in  this  case  a  point 

scatterer)  varies  as  a  =  vlx^+d^  where  a  is  the  distance  to  the 
target,  d  is  the  depth  of  the  target  and  x  is  the  horizontal  dis¬ 
placement  from  the  position  immediately  over  the  target.  At  each 
measurement  point  off-axis  from  the  target,  its  apparent  depth 
is  i.  It  is  at  a  position  corresponding  to  i  in  each  radar  return 
that  the  target  echo  will  appear.  If  these  distances.  Si,  are 
"straightened  out"  as  shown  in  the  lower  portion  of  Figure  32  it 
can  be  seen  that  the  position  of  the  target  echo  traces  out  a  hyper¬ 
bola.  In  general  then,  the  type  of  signature  that  we  are  looking 
for  in  a  cross-sectional  view  is  a  hyperbolic  curve. 

To  verify  this  concept,  waveforms  calculated  by  Davis  [4]  were 
grey  level  plotted.  Each  waveform  was  derived  from  a  moment-method 
calculation  of  the  backscatter  from  a  circular  cylinder  of  e  =1 
immersed  in  a  homogeneous  medium  of  e  =9  and  a=0.01.  The  depth  to 
the  center  of  the  cylinder  was  2.745  m.  Nine  waveforms  were  cal¬ 
culated  using  horizontal  offsets  from  the  tunnel  ranging  from  0  to 
2.44  m  (i.e.,  0.305  m  spacing  between  waveforms).  Figure  33  shows  a 
typical  waveform  from  this  group. 

The  resulting  grey  level  plot  is  shown  in  Figure  34.  To 
generate  this  plot,  waveforms  offset  to  one  side  were  duplicated 
on  the  other  side  of  the  cylinder  to  simulate  a  traverse  over  a 
tunnel.  As  expected,  a  hyperbolic  shape  appears,  accompanied  by 
a  train  of  echoes. 

To  check  the  possible  distortion  of  the  hyperbolic  signature 
by  sparse  sampling  (i.e.,  too  few  antenna  positions)  the  above 
data  was  re-processed  using  fewer  waveforms.  Figure  35  shows  the 
result  of  using  only  5  of  the  original  9  waveforms  and  Figure  36 
uses  only  3.  Note  in  particular  the  increase  in  background  level 
in  these  plots  and  the  widening  of  the  hyperbola. 

Returning  to  experimental  data,  Figure  37  shows  a  cross- 
sectional  view  derived  from  one  linear  traverse  of  the  Curtis 
School  data.  This  data  set  consisted  of  a  group  of  31  waveforms 
taken  along  the  H  coordinate  of  Figure  18.  Measurements  were  taken 
at  column  positions  -5  to  25.  The  time  span  of  each  waveform  was 
200  ns. 

Note  that  traces  of  hyperbolic  arcs  are  apparent  in  this 
plot,  mostly  centered  about  the  17  position  on  the  horizontal  axis. 
The  vertical  repetition  of  the  arcs  is  a  possible  indication  of 
multiple  signal  reflections  inside  the  tunnel,  as  discussed  earlier. 
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i  =  ./x*  +  d2 

Figure  32.  Derivation  of  Target  Signature  in  Cross  Section. 


54 


1— »  *■» lit VA *J f  **' 


l  l.l-tJVIH  K.  n,  tl. 


OFFSET 

K-0.305m 

8765432  101  2  3  4  5  6  7  8 


30.8  ns 


- 2.2.7 r», 


►•*»***►*► 

■  **►*•’ 

» •  *e*k.fb>  ■ 


0.5  m 


ffifmv  ©L« 

Figure  34.  Moment  Method  Datas  Cross  Section. 
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Folded  Dipole  No  Clipping  Fold  First  Using  5  Waveforms 
Figure  36.  Moment  Method  Data  Cross  Section. 
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Figure  37.  Curtis  School  Site  Cross  Section. 


The  dark  line  running  horizontally  across  the  plot  at  the 
top  (labeled  "initial  coupling")  is  the  coupling  signal  between 
the  transmit  and  receive  arms  of  the  antenna  and  may  be  used  as 
a  time  reference  to  locate  the  surface  of  the  ground.  The  other 
horizontal  dark  lines  are  caused  by  reflections  from  the  balun 
transformer  in  the  antenna.  That  is,  part  of  the  transmitted  pulse 
is  reflected  at  the  antenna  terminals  and  travels  back  to  be  reflected 
again  at  the  balun.  This  is  in  essence  a  second  transmitted  pulse. 
These  multiple  balun  reflections  then  result  in  the  multiple  dark 
lines  seen  in  the  plot.  The  presence  of  hyperbolic  arcs  at  a 
shallow  depth  could  again  be  caused  by  the  disturbance  of  drainage 
by  the  tunnel . 

Another  site  from  which  data  was  processed  is  the  Hazel  A 

mine  in  Gold  Hill,  Colorado.  This  is  a  former  gold  mine  dug  straight 

into  the  side  of  a  mountain.  Measurements  were  taken  in  a  large 

grid  on  this  slope  according  to  the  scheme  shown  in  Figure  38. 

Measurements  were  taken  first  using  a  6  ns  pulser  and  a 
Terrascan  antenna  redesigned  for  deep  probing.  As  stated  in  Chapter 
I,  the  modification  consisted  chiefly  of  moving  the  balun  transformer 
closer  to  the  antenna  terminals  to  obtain  a  selected  range  window 
free  from  balun  reflections.  Later  measurements  were  taken  using 
an  antenna  designed  specifically  for  the  type  of  target  found  at 
Gold  Hill.  Transmission  experiments  between  the  tunnel  and  the 
surface  indicated  a  relative  permittivity  of  approximately  6  for 
this  site. 


First,  a  plan  v'jw  of  the  area  was  plotted.  The  depth  of 
the  tunnel  throughout  this  map  was  between  2.75  and  3.32  m  (see 
Figure  38).  For  a  midrange  value  of  3.02  m  we  calculate  the  position 
of  the  tunnel  echo  in  the  waveforms  to  be  about  (3.02  m)/(0.0612 
m/ns)  49  ns  (not  taking  into  account  the  lateral  distance  from 
the  tunnel).  The  tunnel  depth  span  of  3.32-2.75  =  0.57  m  translates 
to  a  time  span  of  9.3  ns. 

Thus,  we  choose  a  time  window  for  RMS  value  calculation  of 
approximately  20  ns  centered  at  the  49  ns  point  on  the  100  ns  wave¬ 
forms.  This  corresponds  to  a  depth  range  of  2.41  to  3.59  m.  Figure 
39  shows  the  plan  view. 

Since  contrast  is  very  low  in  this  plot,  clipping  (see 
Chapter  III)  is  used  to  raise  the  contour  out  of  the  background 
(Figures  40  and  41).  We  expect  to  see  a  dark  form  running  from 
left  to  right  centered  about  the  10  coordinate  but  this  is  only 
partially  visible.  It  is  puzzling  that  a  similar  shape  occurs 
centered  at  the  5  coordinate  also.  Compounding  this  problem,  we 
note  that  in  Figure  42  which  shows  a  plan  view  at  a  greater  depth 
(59-78.2  ns  or  approximately  3.61-4.78  m)  that  a  similar  large 
dark  area  appears.  This  is  perhaps  a  multiple  reflection  from 
the  tunnel. 
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Figure  38.  Geometry  of  Gold  Hill  Map  1  Site. 
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Figure  40.  Gold  Hill  Map  1  Plan  View 
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Figure  42.  Gold  Hill  Map  1  Plan  View 
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We  observe  at  this  time  that  the  data  taken  with  the  Terrascan 
antenna  has  extremely  high  clutter  levels.  It  is  anticipated  that 
further  processing  of  this  data  will  reduce  the  clutter  level. 

Later  results  using  the  long  box  antenna  will  show  substantial 
improvement  in  this  area. 

Following  a  different  approach,  we  next  plot  a  cross-sectional 
view  using  only  the  waveforms  in  a  straight-line  traverse  directly 
over  the  tunnel  (i.e.,  waveforms  with  coordinate  10  in  Figure  38). 

This  is  shown  in  Figure  43.  The  expected  form  of  this  plot  would 
be  a  dark  area  spanning  the  plot  from  left  to  right  (i.e.,  along 
the  tunnel)  and  perhaps  dropping  towards  the  right  because  of  the 
greater  depth  of  the  tunnel  as  we  travel  uphill.  This  type  of 
signature  is  apparent  in  two  places  in  the  figure.  The  first  position 
is  almost  exactly  where  the  tunnel  location  was  calculated  and 
the  second  somewhat  higher  than  the  roof  of  the  tunnel.  This  implies 
that  discontinuities  caused  by  water  flow  and  fracturing  during 
tunnel  construction  are  the  source  of  such  signals. 

The  next  step  was  to  again  plot  cross-sections  of  the  data 
but  using  lateral  traverses  across  the  tunnel  instead  of  along 
it.  This  form  of  plot  would  possibly  exhibit  the  hyperbolic  form 
discussed  earlier.  Figures  44-50  show  this  type  of  plot  for  traverses 
along  the  C,D,E,F,G,H,  and  I  coordinates  (from  0  to  18)  shown  in 
Figure  38. 

On  these  plots,  the  arrow  with  a  "T"  beside  it  indicates 
the  center  of  the  tunnel  and  the  given  tunnel  depth  shows  the  height 
of  the  tunnel  at  its  center.  A  consistent  dark  area  occurs  in 
these  plots  in  the  general  region  of  the  tunnel  (spanning  8-12 
on  the  horizontal  scale).  This  dark  area  is  not  contained  in  other 
regions  of  the  plots  as  the  depth  approaches  that  of  the  tunnel. 

In  general,  this  shape  appears  before  the  expected  time  of  arrival 
of  the  tunnel  echo.  Note  the  continuing  presence  of  shallow  clutter 
in  these  plots. 

These  plots  include  an  enhancement  to  the  GLPR.  In  previous 
plots,  folding  (absolute  value)  was  done  about  the  zero  volt  line 
as  shown  in  Figure  51.  While  in  most  circumstances  this  would 
not  adversely  affect  the  grey  level  plot,  it  can  be  seen  in  Figure 
51  that  if  the  waveform  is  offset  from  the  zero-volt  level,  distortion 
of  the  peaks  can  occur,  as  well  as  the  addition  of  an  annoying 
"background"  grey  level.  To  remedy  this,  the  folding  process  was 
changed  to  fold  instead  about  the  mean  of  the  waveform.  As  shown 
at  the  bottom  of  Figure  51  this  results  in  true  representation 
of  peaks  and  elimination  of  the  background  grey  level. 
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Returning  to  Figures  44-50,  to  aid  in  detection  of  a  possible 
target  a  hyperbolic  curve  for  a  point  scatterer  at  the  appropriate 
depth  is  shown  in  Figure  52.  The  curve  has  been  drawn  to  the  same 
scale  as  the  plots.  This  calculated  curve  should  be  used  with 
care  however,  since  the  actual  target  in  this  case  is  not  a  point 
scatterer  but  a  square  tunnel.  It  can  be  seen,  though,  that  curves 
that  approximate  the  appropriate  hyperbolic  contour  are  present 
in  Figures  44-50.  This  is  indicative  of  the  presence  of  a  tunnel¬ 
like  target. 

One  feature  of  grey  level  plots  (or  any  display  of  underground 
radar  returns)  which  makes  them  more  difficult  to  interpret  is  the 
presence  of  "clutter";  that  is,  targets  in  addition  to  the  desired 
one.  In  underground  radar  some  clutter  may  be  caused  by  surface 
roughness  or  ground  inhomogeneities  or  other  factors.  Since  this 
type  of  signal  should  be  uncorrelated  from  one  measurement  position 
to  the  next,  averaging  of  consecutive  waveforms  in  the  mapping 
traverse  may  improve  the  appearance  of  the  data  by  removing  the 
clutter.  Note  that  this  use  of  the  word  clutter  does  not  include 
actual  subsurface  targets  other  than  the  desired  one  (i.e.,  so 
called  "false  targets").  This  clutter  reduction  method  will  only 
remove  those  signal  components  that  are  uncorrelated  with  traverse 
position. 

This  technique  was  implemented  by  averaging  groups  of  3  adjacent 
waveforms  from  traverses  across  the  hill  and  using  the  resulting 
averaged  waveforms  as  input  to  the  6LPR.  Figures  53-58  show  the 
results  from  applying  the  technique  to  traverses  C,F,  and  I.  Figures 
53,  55  and  57  present  the  unaveraged  data  and  Figures  54,  56  and 
58  present  the  data  averaged  for  clutter  reduction.  Note  the  smoother 
appearance  of  the  averaged  plots  and  the  absence  of  some  of  the 
dark  areas  present  in  the  unaveraged  plots.  Presumably,  the  areas 
"missing"  are  clutter  signals.  Note  also  the  hyperbolic  form  which 
appears  rather  strikingly  in  Figure  54.  This  shape  is  almost  exactly 
that  of  the  calculated  hyperbola  for  that  depth  (see  Figure  52). 

Some  of  the  figures  in  this  group  also  used  another  added 
feature  to  GLPR,  namely  the  ability  to  clip  white  as  well  as  black 
levels  to  improve  readability. 

Comparison  of  waveforms  used  as  input  to  the  GLPR  (Figure 
59)  and  the  resulting  grey  level  plots  (Figures  44,  53,  54)  will 
show  that  the  zero  crossings  of  the  original  data  are  not  readily 
apparent  on  the  grey  level  plots.  The  reason  for  this  can  be  seen 
by  considering  the  sequence  of  operations  used  in  generating  the 
plots.  The  top  of  Figure  60  shows  the  sequence  used  in  the  plots 
that  have  been  presented  thus  far.  The  waveform  is  first  folded 
and  then  interpolated.  The  interpolation  will  tend  to  fit  a  smooth 
curve  to  the  samples  and  as  a  result  the  zero-crossings  of  the 
original  signal  are  "smoothed  out".  If  the  sequence  of  operations 
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Figure  57.  Gold  Hill  Map  1  Traverse  I  Cross  Section. 
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Figure  58.  Gold  Hill  Map  1  Traverse  I  Cross  Section. 


is  changed  to  that  shown  in  the  bottom  of  Figure  60  the  result 
is  different.  This  method  does  not  affect  the  shape  of  the  waveform 
near  the  zero  crossings  so  they  should  be  more  apparent  in  a  grey 
level  plot. 

This  fact  is  readily  verified  by  Figure  61  which  shows  the 
result  of  using  this  FOLD  LAST  method  of  processing  in  the  analytic 
data  for  the  circular  cylinder.  The  hyperbolic  shapes  are  now 
clearly  separated  by  a  white  area  (compare  to  Figure  34).  Figures 
62  and  63  show  also  the  effect  of  more  sparse  sampling,  as  shown 
earlier  in  Figures  35  and  36.  Note  the  reduction  in  background 
level  in  these  FOLD  LAST  plots. 

Experimental  data  processed  using  this  technique  is  shown 
in  Figures  64-73  which  present  the  same  data  as  shown  earlier. 

Figures  64-70  correspond  respectively  to  Figures  44-50  and  Figures 
71-73  present  the  clutter  reduced  data  shown  before  in  Figures 
54,  56  and  58. 

In  observing  these  plots,  it  should  be  noted  that  the  zero- 
crossings  (white  areas)  serve  only  to  delineate  polarity  changes 
in  the  waveforms  and  should  not  be  focused  upon  as  possible  targets. 
Only  the  darkened  areas  represent  strong  signal  returns. 

One  feature  noted  in  the  grey  level  plots  presented  thus 
far  is  the  presence  of  a  large  number  of  echoes  coming  from  scatterers 
closer  to  the  surface  than  the  roof  of  the  tunnel.  This  is  perhaps 
explained  as  the  result  of  reflections  from  locations  in  the  tunnel 
overburden  where  fissures  filled  with  water  or  air  are  present. 

In  certain  cases,  this  clutter  could  conceivably  prevent  detection 
of  the  tunnel . 

One  traverse  was  also  taken  across  this  section  of  Gold  Hill 
using  a  different  antenna.  Tins  was  the  long  box  antenna,  designed 
specifically  for  this  type  of  target.  Data  was  taken  at  0.61  m 
intervals  instead  of  0.305  m  as  before.  The  measurement  positions 
as  identified  using  the  coordinates  in  Figure  38  were  C2,  C4,  C6,..., 
Cl 4 ,  C16,  C18.  As  a  preliminary  step  in  evaluation  of  this  data, 
all  256  samples  were  plotted  in  a  cross-sectional  view  using  the 
FOLD  FIRST  technique.  The  result  is  shown  in  Figures  74  and  75. 

These  are  two  distinct  grey  level  plots.  If  an  attempt  were  made 
to  plot  all  256  samples  in  one  graph,  either  the  data  near  the 
bottom  would  be  too  light  or,  after  clipping  to  raise  the  grey 
level  at  the  bottom,  the  top  would  be  too  dark.  These  figures 
then  consist  of  two  halves  of  the  waveform  set,  each  clipped  at 
an  appropriate  level.  The  difficulty  with  this  plotting  method 
is  that  two  different  interpolations  are  being  presented  and  they 
may  distort  the  data  near  the  border  between  the  two  plots.  A 
possible  remedy  to  this  situation  would  be  to  multiply  the  data 
by  a  linear  gain,  increasing  with  depth. 
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Figure  64.  Gold  Hill  Map  1  Traverse  C  Cross  Secti< 
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Figure  65.  Gold  Hill  Map  1  Traverse  D  Cross  Section, 


►  VI 


B 

5 


93 


ttlBBM^WlAWWWIIPWJ^ITO  V  MK.'V  l'  H‘ «.«, w  ha  w  w  ex 


£ 


-*\  |^— 0.305m  ^T 


32.5ns- 


0 

i  i 


2 

i 


i 


i 


6 

i 


i 


I 


10 

I 


I 


12 

i 


14 


16 


18 


i 


i 


i 


i 


.  ......  .*•••«••«»><  i 

•  v  ... 

,cfrPf.T*tkHMIUIMoTf 

*»»»>* 


...  -1.98m 


II6 


* ,v.t« • .. 

M<v«n 

-  '  *pr»rr r"#  * '  *  * 


i 

_ ,T:\  < . 

♦*  -  ■'*  *"  *»IW*M*»**M*H»* 

>««•»••  4  ♦**■»•. *t »p.  » 

...  -  . . . 

4*.  •»<«»•»»•**  a  -  PAQ 

. UiCDlTl  F  VW 

......  t .rTiim(miBtmiiiiii«i<itmmMiimif(iiJimii 

...  .  .«IIBIBIIBlBBaillll|||ll*<flllBIBBIIBBBBB(BBMB1lll< 

.  #l«fHBIIIBHIII<l»««MI>M»<BffBtHtlMMfMfM|r 

.„<.««pilmHI|BB|l1.<H  lf|BIBBMMBUII»B»MM»Ptrrl#»B|BBB»IBBBMBB«IIB^  * 

•  ..•Mpr*ltUT|l1«ll*BIBtf.«r  I  r  * . .  .^BlllMBBIBBBBBBBBH 

4  .»« »  ,4MBIIBBBBBBIBB|lBlr  »  —  trVB||«l»l«M*IMrffp. 

•  '"Ctrrrr  1  pppBBIMBBf..  •-  .  •.tlllUBMUBBot.'.. 

Iiprmr't"  »«***•  »«*••••••**  -  ♦  •  r**'r'  ••lllMMMBrl  ‘  »*  _ 

K.MO.  ..**p**»**lr  i*  — •  .  .*»r»r*a 

4«i.pl|i(rti4*  »i*«iuuai»*(  *‘r*>t* 

,4.|.nvr.'.«i'* H4M<* '■.n,«mr'r4.-..  .>t«ii«*uiin'»rp  .mivaMATj-v 
**•»  .leitivp'i*.  .«4«ki<4  ... . r4gaaiiiBBBiut;>4 

4<PB»»#«f|kt«nt*4  ••«*!•  .MMlBMIBalr*  4»»p»Bl»MB*BBIB»B#.»P.*»M»tll»*4< 

•»i4fB»IH*tt««kt>aet4  ..«*  *lWlBBIBI»HI»PWr»fBU|BBBMIBMMBt»^f«MM»«,»*' 
.•uklBHlMlMHt«Ut«  .rIHBBIB»IIB«M*.HIt«fBBIBHIBBBH«#l«»BlBIMM 

..•«,rHflBBMBlBMI*«>t-  >(•  «.MBBB»BIMBMMmi*H.I*BIBBB«B|BBH»ltM|BBB4B.T 

•4*>  .  4  ,  f  ►««•»*•»  ***gl>*«tf  •••411*4ptMflBI|llm*l««*'<>  •  • 

■  4  4> t, p»l|f B(||Bll|ll>l«1p'< •’  .4,k»ll(#M1«kl»t1li»4fcM<*4«**-'*l»BIW|MMW*»»w»f* 

. .rgfMMlIlMlIBBlp1*  44iittfB*«l*ttCtl*lt<l  »ar*  X.. 

. . . . . . .  .»»*»».  •ir4T(irM(|«e'<**>_ 

t*.. af. ■-4*«st«*y*.  .*tp«Tr»* 


■M 


8 


-  w 


(•••saBi. 

'rflft" 

*p*» r*er ' • 
r«T*»f  >»»• 
,.p*.«ar'  ' • - 

•>  * 

»• 

•  •  t  »♦ 

. . 


xVMtfetfetcnMB* 


iiT*r»wMlHHB»nkpP»'«erei-fTf* 


'  .IPgcMP1*. 

*  •  - 

P4i^|»»»’»' *« 
1  *••*?**•»•*  *  " 
P  l%***a»»*>'  *  ’ 
r»>****m1*~  • 


4  *  (*»»MI»»P4»*PVPPr«4 

•  >W(tMIIBlBirf1UtpI‘ 

.4«rMrpi»lTIBnp>»p'4 

•  •• 

4*4-.  »♦ 4  » 4  t f f 4  *  ». •« 


.  •  f  4. 


•IIM 

>.*•44 


4<»1**Mf1!fH*lffllB<ll«BIIBMIBB^**lt**M*<fatr4 

4 1  IM***‘P*a^t  4  »MMBB>>«4M^#**'B<p*  * 

•  ftlBBI1**tMI**«i  pglBf|IMP»*4iif*«'4p  t  4 

.ti.Tttf»TTtTTfciiTT*fHTf  f * k T** •*•*»•***«**••« T-.' ; 

4* rp*««*Hrf crrrrrtcc  ‘  b * '  *  .4p|»iwa»M»«*»B»*p  v»4 

4<}*t<ininiMi|Pr<i4 *«».4»  'p4tfTI|«allt|4<*  <  •  • 


vf.r 


p»H»€  I  *• 

4Bg*r<t>ia><M4aki 
»  *»  M»»* 


kyv« wvf v*4~...  '  -*4.ie«»rH  4.  • 

. . .  •  ...» 

. . .  tll»lft4f  f*«. 

k««V»> 4444M I  k«l«4*llf.4(tt 1 *»••♦• 

•1 II »»»f «• *4  444  1 If  r*«|tltSS«a*ccl J »»•♦* 
'*  *#**e  f*”'***  «“«♦♦**•  se*»TTttff*fcep**Trnr4 » 

,  •  ••*(*•  mi«»i4»«4*4*'irP4*t**4M4HP*'  •*• 


:L_TUNNEL 
•  DEPTH 


ppfpf 


10ns  r 


•V '»»♦«*> »>•!*• 


67.3  ns  — 


\  ri-ynxiH 


—  4.11m 


Modified  Terrascan  Antenna 
NXDIV=5  NYDIV=0 
L I  X= 3  LI  Y= 1 


Clipped  at  1.75 
Fold  Last 


-~\  )■*— 0.305m  ^  T 


32.5ns- 


4IIIIIIMUMIHHI*' 


-  1.98m 


.•.(IIHIUIIlUinir.llfHtHil  »»***  ♦  . . *►«*-  -  . —  *  .  • 

•immimniiiiifMiMTR.  r".  . •  •  •  -♦-»  ■  •err*'1 

r\ .  .»tta***Ti-r ■ • 

.itllttlllllllltlltl*.  “•pr*.”- 

fceecePTP’r’'  *?♦«▼•*»»  -• 


. . .  .  >ippi«ifii*m<i.r  . 

-Prr....-*  ....  .. 

»•**»•**>.-  —  '■* 


0.25m  FOR 


...»»»«»  «y  v  v  VMS  tr«* 


r«*.  wv/-  . ***»»  f 

u  M •{ .{ I  <  t •  r««iaNIHMMIMI(iR4r 

r-  **fc,T*Tvv \j.  — . i--  .. - M«i4Pni«7';4i«ria*iMii./M'r'r«Mppp»t<m«MM*(«M<Miia*p 

•  *»•-  .»<(>»-  .»«•*—  •«MM|imiMvrp**tl*iP**»atUMIMHIlllllll«>- 

**■■*•••••«* ♦  •*«•>  «!»»>«.  •»*..  ...il,T«*|l|>l||t<k4tmit.|TNII|IIISHIIIMIIHK4. 

**«••«••■*■# tv.-  *tt  «tvt nfcY* - — nrkttiKfniiiniiiiiMHiitfnTamMHMiiiHMiinkr 

•  Iilllllimtfx  ♦*-  MMMtl*  .M»>itriii|IIIIIIHI|»<ill.iill|()IIIUIHIU«llll|«ir. 

*i*»anMaRi»'«*  •  -  -«i.rrc« ‘  •»»*  ■ e*' *•*•••••••••••-  ’••• 


.<««Ttn<rk'.ie4iiiffik,' 


tefTlNlflli*. 


‘*»iun»*  .  —  •••♦•  -«»►•*»»- 

ifTCk1pf»'  **>*i*«»rxrttr<** - »  *Yr«TTT*t  r 


.tFrMa«K*ti»Maaf«>m>^'r»i«tfvtk->T4*« — i»u>nrnmnu*  - 

..v«MPl.r^Mii|iiirmm«TMuifirr* - ^iTryTtT\r»i*«Mif»»tpree»irtrr*aaMMiHmMifnnr.v 

•  •«»»!•  r  at  •*«asiv«Msh***r  «•*»*■*•>«  *--•«**»> i •  ttf*? tvf«iM«tMiiuif ai<  -  • 
**f RipSlfiMNSStr r ma«h«t» ' . 

•  . .  (im  wmn  r*-** - -i.-{v,rT-«?PPs«ftYArlv « ft^kwtarpmi* - 

.«  .  -.<>  (|  )»•..  «.  -  rer«*r  *  * 4  • 


•  er<iu«i>cr«fiTrTio^n«iii.r1ifiioriki«.. 

akt-pppv»<tM«*iN*i«*»ikk>faaaiiaaTt>iti.v' 

''prrtrMtMMaaaaaaaaMtfa^aMfaat*****— 

.M»44»«>r4«a«aMNi(««lfata«aa>Mft.^** '  — 
<  aP««pf»Sf  *Mtf  p* 

----  -kXfMPMaweiwyi^ee.ii-^^mj' 


...... 

—  •4fiiaaa«i* 
•M»p*#aai«*«i* 


"  *  ”  •**»'■**'  •  - - •‘THkt'Y*. 

.....  ....  ..«■!«> 

..•«!  «y*..  . 

♦*»  irctrct*  .  rcrpf*»»»»*a,M»«M  •— 
-•  *«*..**»  ir)>pf(fi  n>rp**»»*»trM«M*ir  . 

.4Ktki>atiy«*t»umrfa«(tirr»t»a»aTf«aa«»r.«r» 

. . .  <i»prt'  •> 


•'..at*.  .pttp 
•  »*tP 


"TUNNEL 
.:  DEPTH 


4  tm  ttm  *♦  •  »«»...■“ 


—  *7v  iiMTi  vx»*V*i*»  »~ 


4.-^; - 1^44»>*  ' 


67.3  ns  - 


• - -i*\i  \  ,r4»r.-.i.4i,44.;*—---r 


*«.tu*vyvTunTTTTm.tr 


Modified  Terrascan  Antenna 
NXDIV=5  NYDI V=0 
LIX=3  LIY=1 


-4.11m 


Cl ipped  at  1.67 
Fold  Last 


Figure  67.  Gold  Hill  Map  1  Traverse  F  Cross  Section. 
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Figure  68.  Gold  Mill  Map  1  Traverse  G  Cross  Section. 
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Figure  69.  Gold  Hill  Map  1  Traverse  H  Cross  Section. 
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Figure  70.  Gold  Hill  Map  1  Traverse  I  Cross  Section. 
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Figure  71.  Gold  Hill  Map  1  Traverse  C  Cross  Section. 
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Figure  72.  Gold  Hill  Map  1  Traverse  F  Cross  Section. 
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Figure  73.  Gold  Hill  Map  1  Traverse  I  Cross  Section 


Figures  76  and  77  show  the  data  plotted  in  a  range  window 
of  approximately  2-4  m  (32.5-67.3  ns).  Figure  76  is  a  normal  grey 
level  plot  and  Figure  77  is  a  plot  using  the  3  waveform  averaging 
technique  for  clutter  reduction  discussed  above.  Figure  78  shows 
the  shape  of  hyperbola  to  be  expected. 

Figure  76  exhibits  this  expected  shape  at  a  position  less 
than  10  ns  above  the  roof  of  the  tunnel.  We  may  also  compare  Figure 
76  with  Figure  53  which  presents  data  from  the  same  traverse  across 
the  hill  but  obtained  with  the  Terrascan  antenna.  Because  of  the 
different  antenna  structures,  the  clutter  should  appear  differently 
in  each  plot;  and  the  figures  show  this,  but  note  the  darkened 
area  that  appears  in  both  plots  near  the  tunnel  position. 

The  spreading  out  of  the  hyperbolic  shape  in  the  plot  using 
the  long  box  antenna  can  be  explained  in  two  ways.  First,  the 
spreading  could  be  due  to  the  lower  frequency  response  of  the  antenna. 
Secondly,  the  traverse  taken  with  the  box  antenna  had  antenna 
positions  separated  by  0.61  m,  or  twice  the  spacing  of  the  Terrascan 
version.  As  was  shown  earlier  in  Figures  35  and  36,  more  sparse 
mapping  results  in  broader  hyperbolic  forms  in  the  grey  level  plot. 

Plots  were  also  made  using  the  FOLD  LAST  method.  Figure 
79  shows  the  unaveraged  data  and  Figure  80  shows  the  averaged  data. 
Note  once  again  the  smoothness  of  the  clutter  reduced  plot  and 
its  tendency  to  single  out  one  or  two  hyperbolic  forms  as  candidates 
for  consideration  as  targets.  Observe  that  the  averaging  process 
tends  to  remove  the  rapid  variations  in  null  positions  which  are 
also  clutter-related. 

Other  data  from  the  Gold  Hill  site  was  also  subjected  to 
grey  level  processing.  Two  diagonal  traverses  were  made  using 
the  long  box  antenna  at  a  location  further  uphill  than  that  of 
Figure  38.  Figure  81  shows  the  layout  of  the  measurements  taken 
at  this  site.  These  antenna  positions  were  the  only  ones  available 
in  this  region  because  of  the  rough  ground  and  presence  of  many 
obstacles. 

We  will  first  discuss  the  traverse  from  L-2  to  S12.  The 
complete  256  point  waveforms  are  shown  plotted  in  grey  scale  in 
Figures  82  and  83.  They  are  two  separate  plots,  as  discussed  under 
Figures  74  and  75  above. 

After  selection  of  an  approximate  depth  window  of  4.6-8.84  m 
(or  75-145  ns)  for  the  tunnel  depth  of  6.1  m  in  this  area,  a  variety 
of  processing  techniques  were  used  on  the  data.  Figure  84  shows 
the  basic  plot,  clipped  to  highlight  the  areas  of  interest.  The 
calculated  hyperbolic  shape  for  this  map  is  shown  in  Figure  85. 

Figure  86  shows  the  same  data  after  clutter  reduction  by  averaging. 

Note  in  particular  the  noise  evident  in  the  lower  portion  of  these 
plots  and  also  in  the  typical  waveform  from  this  traverse  shown 
at  the  top  left  of  Figure  87. 
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Figure  77.  Gold  Hill  Map  1  Traverse  C  Cross  Section. 
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Figure  82.  Gold  Hill  Map  2  Traverse  L-2  to 
SI 2  Cross  Section. 
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Figure  83.  Gold  Hill  Map  2  Traverse  L-2  to 
SI 2  Cross  Section. 


112 


w 


iS 


< 


0.683m 


75.0  ns- 


L-2  MO  N2 

J . l . It, 


T1 


145  ns- 


. . 


|T 

RIO  SI2 
*1*  !  "I  l-4.59m 


fe^eMfee^re*****^********#*******"**^****?  »**c»-p**"*  *•••**•■»*••«■•  0.5m  FOR 


€r  =  6 


. . . . 

i(i,m«t«a«ati«ff<e(<i<>««»**»**.*«»*»M*««»»»««»«iiia«creiiiitrMMaaaHa»aivtiiii(tiiMt 
v»aM>i|aaaaaaMw»i««n«<*>,>>'in«***»*****,>***«**,>|»**<,tcr*|l,«»Mnanaa<  M«»«errmi 

»*«PCa»t>«l  Cp«  %%» i*  »*  »i  x  x«Mrrea"**v*»«a«************t»»r  *f  A*. 

“retrer»rjer«M»n«M.#*kn«ir»»*ift«M«M*arraaa»«a«ir»*E,«« 

l.««4*<*»»»M«<iimmii«firrpr*a»<i«***erMi)><**^**«ii?fee*''*a»«ai*iaaiaaa«aiiiHfkitii 


..***u««itt('tc»mtfekr«*a«»cr’rFFrer#«iaaaa*w»aa«afMiMaafa«na»»»»Mt»»»(»- 


l.«««»«t>ku>m(mart(trtu<it>*,«»..«Mi»r*»i«M«aiiaiiiaaiaaaiiuaiaiaaiaaaaiika*ti«t<> 

...  ...««»>iiiitMktrttet!kkit*aaaaiaaaiiaaaaaiaaaiMaiaa«aaiaiaaniiaiianafM«i«t 

. ^^inttirrro.HctMaMiiaaMaiaaaiataiaiaiiaiiiiaiiiiiMaaaamvMiKh 

.  . . 

...»»MiuivM*«.«**M«4«iitR4Tiii«(aaaiaaaiaiiaaaaaaaaaiaaiiaaafaai*aaMa>aai<iiiMWM« 

..«**.x  mix. .......  •••iMiiiiiaaiaiaaaaaaaaiaaMiiiaaaaiiiiiiiaaaaaaiiiiamuaaaaa* 

. . . 

,,,»iiimi>ivi>miiiii>irc*H*aauiaaaaaaciaMiaiaiaaa>iagaiaNaiiaiaaaaiaaaa<aaaaiiaiaaa 

i«.*4*»t.««.«*M<««t«,(iH«MaMiMMaHiiaiiaiiiaiaaaaaaiaaaaaaaaaiiiiaiaaMauaaauaiaiaa 

nriaai»«iiirm*»>*ifaaaaaaa#MaMNaaaaf>fr*rfett»a«)iM»»««»«*«ns  Ttt»»i'**********« 

,.P*  «**<»  «■  ■  *>»•  t « ^pc.«  •  t  ...... ..-mi*..  ♦  «»  *  ■  ■n«»veeerc*1*a** 

iMfenttiPkerMHOMiiMrr»<*u«»**.>»«Miiiiertii«tteei*  »*x •♦♦♦♦•♦•♦•»**» it t  ittwunpit  i 

luUXMHitttOnt^Ha^eFMrFtiTM'aanrtvetc^tectMiiM**— ..«**«nii>utttH*tI  r 

>*a««»»rcecc»*  «•*>***•••♦«•  A*  <«  ter 

•  a  »♦♦♦.—.  ....« i  >  <»«itMaiMUaiiHHlaiiNNa»aiaa«f  kin  mi  »•*•••• 

....  >Maaatwaaiaa«aiilni>n>aaaaa»aan«»ti»  . . . 

>t1i1««....—....t*M(itppiifMlpiaiiaiaaii<a<iiaiNaaHffaifMp<itu<kiH<i 
.  ,Y < " i <* c«r r««t»#i##i#»*Miiaaf j ,»f • ,  ><>ui***»»»*  u 

luinmmMitppotiiitaaanaaiiiiaaMaaaaaaaaiMaaiHiiHMxatMiititiiiteriiitumiu  - 

. . .  *••»?*«  •••■•ccerce  **»*»»»  x  »•* 

. •.>4t4re.»mo»MMfM»aaaaa»aaaiiaMiaaiaaaaa'ia«i»M«**»M»mii»t»»»fert>* 

..*»*«>i«rrmia«i«<Mut»«Miii(aaaaiaaaaHwaa<MMaM*MaHMMMMiiiMife» 

k««rr»**araM«M*»»«aaailHNni«a«at«aiia*rcrrtcep«»a«*afafHri|MiM* 

••••*«•«.  Tf  ft  T I *««*•*•  •••••••••fTTf*. 

. .  • 

. . . . .  - . ......... . . .  ' 

. . . . . . -  - *««*•*  -.p.*rrr rtii rrr/irr:>«a.n*»»*-«.r 

. . «•  •*"»*••  - - ............... - iti 

• 

—  ..•««**.  «« t  •  i  tt  i  .>*«**.•.... — . ...  > 

. . .  reccreteMerprrt***  iti  M*»««»«t««attm«( 

f  *t  •►■tPtPtPt  •*••••*»»►  •■■r»pr*-tee**-pprrrctpftepcrrpr-;ec 


TUNNEL 
DEPTH 
AT  RIO 


mTf|M»tfarr* 


***’  p«x*  •«•••■*•  *«•»*•*.  'tnUMfjijt.mt.npjfajjM 


i»«»««iiH*aaaaaaHaaiaaaaiaiMmM«M«tito»»p!iuu4i  i  j  ••  •  . . . 

•  **4««4»*4  ■ 1  *  •  •  tce**«*<tatiiatiiiitiH*fH«aaitaaMMtfa»)i>4  •  ■  ♦  ♦»«* 

•  *  »  1  .*».....  .♦>  1  »**tl«*liai||H|PMlttt»m«IMM4a 
»*«>«iii»«rrrM*««i|TmntilTimYMtafifittmfti<M.rcrr  ittttttittrrcai* 

. it  tax o  *»> 1 t'crr*" 

uieetPai»M*«iMMMM«atiao»Nanaa««H*aaiMauMiMiiiiiiiiaauisiftii 

. . . .  i.«t«k(*.|l|||ttlt|aia 


—8.84  m 


Long  Box  Antenna 
NXDIV=1 2  NYDIV=0 
L I X= 3  LIY=1 


Clipped  at  1 .5 
Fold  First 


Figure  84.  Gold  Hill  Map  2  Traverse  L-2  to 
S12  Cross  Section. 
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Figure  86.  Gold  Hill  Map  2  Traverse  L-2  to 
S12  Cross  Section. 
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Determination  of  Best  Trapezoidal  Low-Pass  Filter 
for  Gold  Hill  Map  2  Data. 


Low-pass  filtering  was  chosen  as  the  method  to  remove  this 
noise.  That  is,  the  waveforms  used  as  input  to  GLPR  would  be  filtered. 
Figure  87  shows  the  effect  of  several  successively  narrower  low  pass 
filters  on  one  waveform  from  this  traverse.  A  trapezoidal  filter 
was  finally  chosen  with  lower  corner  frequency  at  100  MHz  and  upper 
corner  frequency  at  140  MHz.  This  filter  removed  most  of  the  noise 
without  adversely  affecting  the  waveshape.  Figure  88  then  shows 
the  grey  scale  plot  of  the  filtered  data.  Figure  89  results  from 
first  filtering  the  data  and  then  averaging  it  for  clutter  reduction. 
This  clearly  shows  the  expropriate  darkened  area  in  the  region 
where  the  tunnel  response  is  expected.  It  also  contains  the  response 
from  both  the  roof  and  floor  of  the  tunnel. 

The  data  from  this  traverse  was  also  processed  using  the 
FOLD  LAST  method.  The  same  set  of  processing  techniques  was  used 
as  before.  Figures  90-93  show  the  results  of:  1)  just  clipping, 

2)  3-waveform  averaging,  3)  filtering  and  4)  filtering  and  then 
averaging,  respectively.  Note  that  in  this  group  of  plots  there 
are  two  hyperbolic  forms  at  approximately  the  right  depth  and  that 
they  have  the  same  general  shape  as  the  calculated  curve  (Figure  85). 
This  could  perhaps  be  explained  as  reflections  from  both  the  roof 
and  floor  of  the  tunnel. 

The  next  group  of  figures  presents  the  data  from  the  other 
diagonal  traverse  shown  in  Figure  81,  namely  that  from  0-2  to  Xl6. 
Figures  94  and  95  are  two  halves  of  the  complete  waveform  set  with 
Figure  94  spanning  0-100  ns  and  Figure  95  100-200  ns.  It  can  readily 
be  seen  that  noise  was  a  problem  on  these  measurements  also  by 
observing  Figure  95. 

The  depth  window  chosen  for  the  last  traverse  (75-145  ns) 
was  also  used  here.  Figure  96  shows  the  basic  grey  level  plot 
while  Figure  85  is  the  expected  hyperbolic  shape.  Figure  97  is 
the  clutter-reduced  version. 

The  same  low  pass  filtering  method  (with  the  same  filter 
corners)  was  used  also.  Figure  98  is  the  unaveraged  but  filtered 
data  and  Figure  99  is  the  data  after  first  filtering  then  averaging. 

Versions  of  these  plots  using  Ihe  FOLD  LAST  method  were  also 
made.  Figures  100-103  are  the  0-2  traverse  after  1)  normal  grey 
level  plotting,  2)  3-waveform  averaging,  3)  low-pass  filtering 
and  4)  filtering  then  averaging,  respectively. 

Another  useful  processing  technique  was  also  tried  out  on 
this  data.  As  shown  in  Figure  104,  the  normal  use  of  the  folding 
function  (i.e.,  absolute  value)  in  grey  level  plotting,  using  either 
the  FOLD  FIRST  or  FOLD  LAST  methods,  causes  positive  and  negative 
peaks  to  be  mixed. 
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Figure  88.  Gold  Hill  Map  2  Traverse  L-2  to  SI 2  Cross  Section. 
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Figure  89.  Gold  Hill  Map  2  Traverse  L-2  to  SI 2 
Cross  Section. 


USB 


&$ 


I 

m 


pi 


V.V.V. 


m 

Wi. 

S®&: 

MW 


<  J v.% 
w 


0.683m 


L- 

75.0  nsJ. 


i¥fpll.  .»  . . . 

. . r>Mrrr(it.''>i>V)*^M1ffHIIMIIIIiaitlli«tm«ii<ircte 

•***■•*'  •  -»  ••{.»>••  If  Tt  •!  > 

«•«>•«•*  ..  -•  (,pr»l 

»  .  ......TfffiMMMM-iMirMIU'milllllHMIMIIimM. 

»■***«•'♦  -•»  »,,pr»4i**TTfkT*4«>44#*4SMfSSS4i*lM'4*»-»T*«SS*S*M»STlp4 

#.<.<>•  -*•  career »• . «» *vtrrr**«»*«‘*r' 


:;:...:;r.;.::;:;..-4.59m 


•  -•«  •  :»TTMT«»CC  '  •- 

.  .<r  *•*■•■■■•»»•**»?’  * 
tieMrlllMlPftlt'i 


0.5m  FOR 
«r=6 


. . .  „  w*’-r _  *■ 

,,  . .  . . 

,....rrP  . . . .  . 

..  .,..rr  i  »  •  »  X*  _ _  - 

. — — 

ii  (•<•»—•  ixu>s»t«*«Kcccfckt 

feifti«c)UiXi*i _ -  — -  - 

.4,^«  . . . . 

. . . 

. . . 


i  >  »»  «**♦«»  <*r>rccrfcci-rtc'*  ••«*<!  t  »ff  l* ■  •♦« 

II  kM1|lilk»»*««  •.♦llifCimlktMklOliMHUlMkHlC*  *  '* 

....<..cccct;ccetccl..jXi,  ..  •.>Xltck'nilHlliniMHMIUlllff«!'rl* . .  - 

..iinimMitiiu.***  ..i  »r 1 N  * 

-«**«  ti  i  «i«s  /»«*♦- 

.  . .  .Mlif.iiiMiUiliaillfIMHiHMilillliiilHMtnetxif.... - 

...■llltltliT**.  ..■U.llMlIMlIUMdUHMfUIMIIIimilUtMMIliirFt'ii''***** 

. . —  *♦**- 

..i.tiVMllHIIlHIIIHIHMMHMlMMOIttMltlfHvKMPIIfikCtP11'4** 

. . . 

. 

—  . . . 

.......  «.  ***? 

.  . . . .  n»m“ 

. .  •«  ...*♦■•>#»•»**«*  ererer  ■ 

4rlkHt|tt'*"v  ♦♦•••  ••44Mi,,'I»P*',MI>;*»*«  «-««•«  «k  1 t  ll  **»*4»i|  »».* 

i*i  iptf.cr'**«**a*  ixx •♦•••• 

• .  i*»i  • ......  . in  *  •  tktxiitrm  itiitnt  »i  . . 

•  »«w  Pt*  V  . . . 

. .  •4i>kl.tMI»tlllll|lll«*'lfM>»ltM|t<llmii.44 

....  •.iMrr»MliMt«HIMNUWNM»«li|t.*kf  X»  *  «*♦ 

......  i  <  i  <  M«pmn  •  »♦♦»♦♦♦♦♦♦♦♦♦♦♦♦•— — 

.......  .EeilPMWHMIlHIIMMMIlHHXWIT^IPCEtElllUM « IB »■»> - 

.  •n.(,tir«4».»»<M*M*».»»»»*'l4rr  . .  •  •rf..*..»4.»M4frf. . "*** 

.  ....  . . . 

4  »•  *  »*  4  »*«»»•  •  *«• .  »!•••■  .....MfcM.M.ik.Mi  •** 

................  . . . . 

***..,,  ...... _ _  .......  ....  •  •  .h».a. «*•*•«». *4  -- 


TUNNEL 
DEPTH 
AT  RIO 


10  ns 


•cteetfett 


perrrt»»*ff  •  'tetitrretrprf  M  f  i 


145  ns— 


. . 

...  »«••»•>  » •»»iMUtereectrr*in»i4‘** 

..  _ ......  ..(,)t{rp.M>«»»Mia*kVyi|MM«fMf1>r>n4.  *I»  »  k  » 

.  . . . 

... ..  »k »,  i4|t<krtvtct*k***k*p(fpffr>****4*l|,kk*p|lk*** 

...,>iirrrr..«.rrrrciittiiPrrp«'PP,'*1'**ll|rr,,>»>  -•ine»t»s«ssss»*s*si««*acc**»4  +  * 

. . - •kUlilJIfftEfllllllltteeliult*. 

..  ...  4  > .  1 1  Mfkterrrtkpr « 

....  . . . . ft <rrrtw*eer*' ».»•••  . .  B8  S4fD 


Long  Box  Antenna 
NXDIV=12  NYDIV=0 
LIX=3  LI Y=1 


Clipped  at  1.78 
Fold  Last 


Figure  90.  Gold  Hill  Map  2  Traverse  L-2  to  512 
Cross  Section. 
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Figure  91.  Gold  Mill  Map  2  Traverse  L-2  to  S12 
Cross  Section. 
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Figure  92, 


Gold  Hill  Map  2  Traverse  L-2  to  SI 2 
Cross  Section. 
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Figure  96.  Gold  Hill  Map  2  Traverse  0-2  to  X16 
Cross  Section. 
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Figure  97.  Gold  Hill  Map  2  Traverse  0-2  to  X 1 6 
Cross  Section. 
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Figure  98.  Gold  Hill  Map  2  Traverse  0-2  to  XI 6 
Section. 


KWMWK 


$ 

$ 

w 


4 

& 


.»!«, 

& 

i*»; 

$ 


ft 

$ 


•X 


Q2 

I 


R4 

I 


S6 

i 


T  8 

I 


uip 


75.0  ft*- 


PO 

I 


V 1 2 


W14 

I 

~  4.59  m 


"T 


;  0.5m  FOR 


10ns  : 


LOW-PASS  FILTERED  AND  AVERAGED  IN  GROUPS  OF  3 

Long  Box  Antenna  Clipped  at  1.14 

NXDI V=1 6  NYDIV=0  Fold  First 

L I X= 3  LI Y= 1 

Figure  99.  Gold  Hill  Map  2  Traverse  0-2  to  X16 
Cross  Section. 
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AVERAGED  IN  GROUPS  OF  3 

Long  Box  Antenna  Clipped  at  1.17 
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Figure  101.  Gold  Hill  Map  2  Traverse  0-2  to  X16 
Cross  Section. 
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Figure  102.  Gold  Hill  Map  2  Traverse  0-2  to  X16 
Cross  Section. 


tmjt  UM** 


7Sj0m- 


WI4 

I 

-  4.59m 


0.5m  FOR 


I::::::::;:::::;:::-:  <r  =  6 


'  TUNNEl. 
■  ■  DEPTH 

AT  U!0 


LOW-PASS  FILTERED  AND  AVERAGED  IN  GROUPS  OF  3 


Long  Box  Antenna 
NXDI V=1 6  NYDIV=0 
L I X= 3  LI Y=1 


Cl ipped  at  1.01 
Fold  Last 


Figure  103:  Gold  Hill  Map  2  Traverse  0-2  to  X 1 6 
Cross  Section. 
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POSITIVE  PEAKS 
ONLY 


If  only  one  polarity  needs  to  be  tracked  to  detect  the  tunnel 
then  plotting  both  polarities  on  one  graph  increases  the  visual 
"confusion"  level.  Plotting  only  one  polarity  of  a  waveform  set 
should  increase  the  ability  to  discern  a  target,  as  is  also  shown 
in  Figure  104. 

This  type  of  plotting  was  done  using  the  0-2  traverse  data 
without  filtering  or  averaging.  The  results  are  shown  in  Figures 
105  and  106.  Figure  105  exhibits  the  positive  portions  only  and 
Figure  106  exhibits  the  negative  portions  only.  These  plots  may 
be  compared  to  Figure  100  in  which  both  polarities  are  plotted. 

Figure  105  would  presumably  show  the  reflection  from  the  top  of 
the  tunnel  and  it  does  appear  to  be  present  in  the  midst  of  much 
clutter.  Figure  106  should  show  the  reflection  from  the  floor 
of  the  tunnel  since  the  reflection  coefficient  at  the  lower  interface 
is  opposite  in  sign. 

Just  prior  to  completion  of  the  present  thesis,  a  new  pro¬ 
cessing  technique  which  appears  to  offer  much  promise  for  future 
mapping  efforts  was  developed.  Using  data  from  the  Gold  Hill  Map 
2  traverses  (specifically  waveforms  T8-X16),  a  set  of  poles  describing 
each  waveform  in  the  Z-transform  domain  was  calculated.  Using 
the  physical  dimensions  of  the  tunnel  to  determine  its  natural 
resonances,  only  poles  near  these  resonances  are  kept  for  further 
processing,  effectively  deconvolving  out  the  unwanted  signal  com¬ 
ponents  from  the  antenna  structure. 

A  new  waveform  is  then  generated  using  only  the  selected 
poles.  This  waveform  should  include  the  tunnel  response  almost 
exclusively.  Some  clutter  may  be  present.  Figure  107  shows  the 
grey  level  plot  obtained  from  this  deconvolved  data.  Note  that 
the  tunnel  height  and  vertical  position  correspond  exactly  to  the 
calculated  values  and  that  the  shape  that  appears  is  very  similar 
to  the  analytic  one  shown  in  Figure  63.  The  only  apparent  error 
is  the  lateral  displacement  of  the  hyperbolic  shape  to  the  left. 

This  is  possibly  explained  by  the  fact  that  position  calculations 
for  the  tunnel  did  not  include  the  side  slope  of  the  hill  on  which 
measurements  were  taken.  As  shown  in  Figure  108,  if  the  traverse 
over  the  tunnel  is  not  horizontal,  the  minimum  depth  to  the  tunnel 
(corresponding  to  the  vertex  of  the  hyperbola)  would  occur  at  a 
position  offset  to  the  left. 
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Figure  105.  Gold  Hill  Map  2  Traverse  0-2  to  XI 6 
Cross  Section. 
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Figure  107.  Gold  Hill  Map  ?  Traverse  T8  to  X16  Cross  Section. 

(An  improved  version  of  Figure  107  is  given  in 
Figure  49  of  Report  4460-11.) 
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CHAPTER  V 
CONCLUSIONS 


From  the  results  presented  earlier,  it  is  clear  that  tunnel 
detection  through  the  use  of  grey  scale  plotting  can  be  accomplished. 
Some  of  the  best  examples  of  this  are  Figures  29,  43,  54  and  107  in 
Chapter  4.  Although  the  detection  capability  appears  to  be  quite 
good  in  these  plots,  the  limitations  of  the  mapping  technique  must  be 
kept  in  mind. 

There  is  a  "learning  period"  in  which  one  must  become  familiar 
with  target  characteristics  as  exhibited  in  the  plots.  In  certain 
situations  it  is  necessary  to  be  familiar  with  the  properties  of  the 
interpolation  method  used  in  order  to  distinguish  between  "real"  tar¬ 
gets  and  those  which  are  created  by  the  interpolation.  Use  of  a 
localized  interpolation  controls  this  phenomenon  quite  well.  Thus, 
in  the  plan  view  type  of  plot,  where  actual  2-dimensional  interpola¬ 
tion  takes  place,  it  is  unsatisfactory  to  use  large  subgrids  for  the 
interpolation  because  of  the  high  probability  that  such  false  targets 
will  appear.  These  are  due  chiefly  to  the  higher-order  polynomials 
involved  in  this  case.  As  has  been  discussed  before,  RMS  window  size 
is  also  a  critical  parameter  in  a  plan  view  plot. 

Preprocessing  of  the  radar  return  before  input  to  GLPR  is  almost 
always  necessary,  as  evidenced  by  the  above  cited  figures.  Techniques 
used  to  generate  these  plots  ranged  in  difficulty  from  simple  DC 
removal  to  deconvolution  using  the  Z-transform. 

One  limitation  of  GLPR  is  its  speed.  The  preprocessing  certainly 
lengthens  the  time  required  to  generate  one  plot,  but  the  real  time 
constraint  is  the  interpolation  speed  and  size  of  map  to  be  generated. 
For  a  YDIM  x  XDIM  sample  array,  processing  a  grey  level  plot  is  of 
computational  complexity 

0(XDIM*YDIM*NXDIV*NYDIV*(LIX3LIY+LIX2LIY)) 


where 


XDIM, YDIM  are  the  dimenisions  of  the  waveform  sample  array. 

NXDIV,NYDIV  are  the  number  of  interpolated  points  required 
between  samples  in  the  x-  and  y-dimensions, 
respectively,  and 
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LIX.LIY  are  the  x-  and  y-dimensions  of  the  subgrid  to  be 
used  for  interpolation,  respectively. 


It  can  be  seen  that  an  increase  in  interpolation  subgrid  size 
causes  cubic  growth  in  computation  time. 

Lastly,  a  fundamental  problem  with  our  particular  plotting 
method  is  the  use  of  a  line  printer  as  the  output  device.  Unequal 
character  dimensions,  uneven  inking,  white  areas  between  characters, 
and  character  skewing  all  detract  from  the  final  plot. 

In  our  use  of  GLPR,  however,  the  above  limitations  have  not  been 
crucial.  Our  interest  is  not  to  provide  full-scale  detection  capa¬ 
bility  with  these  plots  but  merely  to  be  able  to  discern  in  a  simple 
way  our  ability  to  detect  certain  types  of  targets  with  given  radar 
systems.  This  plotting  method  is  also  useful  as  a  measure  of  the 
effectiveness  of  various  signal  processing  techniques. 

Improvements  to  GLPR  are  possible.  Recent  developments  in 
interpolation  using  the  Discrete  Fourier  Transform[9]  suggest  that 
computation  time  could  be  improved.  A  basic  MxN  interpolation  using 
the  Lagrange  method  is  of  order  M^+N^  while  the  Discrete  Fourier 
Transform  method  using  the  FFT  algorithm  is  only  of  order  Mlog2+Nlog2N 
(for  a  rectangularly  separable  interpolation  function,  as  we  have  used). 

Another  improvement  would  be  the  use  of  color  instead  of  grey 
levels  in  plotting.  The  eye  can  distinguish  differences  between 
adjacent  colors  more  readily  than  adjacent  grey  levels  that  are  close 
in  value.  This  is  useful  since  the  reading  of  a  grey  level  plot 
usually  involves  tracking  the  interface  between  two  grey  levels. 
Assignment  of  colors  to  signal  levels  would  be  user-controlled  to  allow 
enhancement  of  areas  of  interest. 
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APPENDIX  A 

DERIVATION  OF  GENERAL  FORM  OF  INTERPOLATION  PROBLEM 

As  stated  in  Chapter  II,  we  are  looking  for  a  polynomial  estimate 
Fp(x)  of  degree  n  to  a  continuous  function  f(x)  of  the  form 

Fn(x)  =  I  akxk  f A1 ) 

We  also  constrain  Fn(x)  so  that 

Fn(xj)  =  f(xj)  j=0,l , • • • ,n  (A2) 

Combining  (Al)  and  (A2)  we  have: 
n  . 

I  akx;  =  f(Xi)  j=0,l • ,n 
k=0  K  J  J 

or  in  matrix  form: 


We  want  a  solution  to  exist  to  (Al)  for  all  x  in  some  interval 
alxlb>  so  we  append  this  equation  to  the  matrix  also,  yielding 


u  m  mm  Kmv.  ro  wwro  w  n»  w  g  u  g  ?  w  w  w  *■  m  g>  ^  i-  m  ^  w  ^ ?■ 


X 

xo 

XI 


X2 

X1 


xn 

xo 

X1 


ao 

Fn(x)' 

al 

f(Xo) 

= 

^ (x7 ) 

an 

f(xn) 

(A3) 


X.-J  a=f 


where 


the  dimensions  of  Xj  are  (n+2)x(n+l) 
the  dimensions  of  a  are  (n+l)xl,  and 
the  dimensions  of  f  are  (n+2)xl 

From  linear  algebra[10],  for  an  MxN  system  of  inhomogeneous 
equations,  X^a=f  such  as  this,  the  condition  for  solvability  is  that 
the  Mx(N+l)  matrix  X2  formed  by  adjoining  column  vector  f  to  Xj  has 


RANK  X2  =  RANK  X-, 
So  we  must  have  that 


( A4) 
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1  xn  x2 
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f(xn)  1  xr 
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(A5) 


vi  y 


la 


We  must  consider  two  cases  for  (A5) 


Case  1: 


j=0 ,1 , *  * • ,n 


Since  the  columns  and  rows  of  Xj  are  independent  (xj  distinct 
and  x/xj)  we  have 

RANK  X7  =  n+1  .  (A6) 

To  find  the  rank  of  Xp  we  use  the  following  theorem  from  linear 
algebra[ll]: 

Let  A  be  an  MxN  matrix.  If  A  has  an  rxr  submatrix  S_ with 
DET  S/0  and  if  every  (r+l)x(r+l)  submatrix  T  of  which  S^  is  a  submatrix 
has  DET  T=0  then  RANK  A=r. 

In  our  case  the  matrix  A  corresponds  to  Xp.  We  wish  to  have 
RANK  Xp=RANK  X7=n+1,  so  we  look  for  submatrix  S  with  DET  S/0  and 
r=n+l . 

We  take  as  submatrix  S  the  matrix 


(,  xf 


1  x^  x^  •••  xn 

n  n  n 

which  is  (n+l)x(n+l).  Since  the  Xj  are  distinct 
DET  S/0 

satisfying  the  first  part  of  the  theorm. 

The  only  (n+2)x(n+2)  submatrix  T  that  we  can  have  is  the  matrix 
Xp  itself.  So,  in  order  to  satisfy  (A4)  we  must  have: 


DET  Xg  =  0 


m 


m 


f(x0)  1  x0  x0 
f(x])  1  x]  x^ 


f(xn)  1  Xn  Xn 


This  is  the  condition  for  solvability  of  (A3)  in  this  case. 


Case  2: 


je{0,l , • • • ,n) 


Since  two  rows  of  the  matrix  _X]  are  equal,  we  have  that 

no.  independent  rows  =  no.  independent  columns  = 
n+1  =  RANK  X]  which  is  the  same  as  Case  1. 

Using  (A8)  we  must  have 


Fn(xj)  1  xj 

|f(x0)  1  xQ 


f(Xj) 


f(xn) 


This  is  satisfied  since  we  have  made  the  constraint  that 
Fn(Xj )=f (xj ) .  Therefore  (A8)  is  the  general  statement  of  the  inter 
polation  problem  presented  in  (Al)  and  (A2). 
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APPENDIX  B 
COMPUTER  PROGRAMS 
PLOTTING  SUBROUTINE 
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